NASA MEMO 5-26- 59L 


NASA MEMO 5-26- 59L 


NASA 

MEMORANDUM 


MODIFIED MATRIX METHOD FOR CALCULATING STEADY-STATE 

SPAN LOADING ON FLEXIBLE WINGS IN SUBSONIC FLIGHT 

By Patrick A. Gainer and William S. Aiken, Jr. 

Langley Research Center 
Langley Field, Va. 


NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


WASHINGTON 

June 1959 





NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


MEMORANDUM 5-26-59L 


MODIFIED MATRIX MEIHOD FOR CALCULATING STEADY-STATE 
SPAN LOADING ON FLEXIBLE WINGS IN SUBSONIC FLIGHT 
By Patrick A. Gainer and William S. Aiken, Jr. 

SUMMARY 


A method is presented for shortening the computations required to 
determine the steady-state span loading on flexible wings in subsonic 
flight. The method makes use of tables of downwash factors to find the 
necessary aerodynamic -influence coefficients for the application of 
lifting-line theory. Explicit matrix equations of equilibrium are con- 
verted into a matrix power series with a finite number of terms by uti- 
lizing certain characteristic properties of matrices. The number of 
terms in the series is determined by a trial-and-error process dependent 
upon the required accuracy of the solution. Spanwise distributions of 
angle of attack, airload, shear, bending moment, and pitching moment 
are readily obtained as functions of qm^ where q denotes the dynamic 

pressure and m R denotes the lift-curve slope of a rigid wing. This 

method is intended primarily to make it practical to solve steady-state 
aeroelastic problems on the ordinary manually operated desk calculators, 
but the method is also readily adaptable to automatic computing equipment. 


INTRODUCTION 


As a result of experience in calculating the span loading on flex- 
ible wings, in which most of the work was done on a manually operated 
desk calculator, it has been found that results are most rapidly obtained 
by combining certain features from several methods. Reference 1 was 
found to be invaluable for its explicitness in defining the relationships 
between aerodynamic loadings and structural restraints in the form of 
matrix equations. Reference 2 furnished tables of downwash factors which 
can be used to eliminate most of the time spent in evaluating the elements 
of the aerodynamic-influence-coefficient matrices. References 3, 4, 
and 5 furnished an approach to the problem in which the overall problem 
was divided into two parts: (l) the determination of the airload produced 

by the various parameters in a particular equation of motion and the solu- 
tion of the equation of motion for the relative amounts of each par ame ter } 
and (2) conversion of the matrix equation of equilibrium into a finite 
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power series which is readily evaluated foi any flight conditions of 
dynamic pressure and Mach number to obtain angle of attack, airload, 
shear, bending moment, and pitching moment resulting from any specified 
initial angle of attack. 

The purpose of the present report is to combine the various features 
of these methods into a procedure which can readily be followed by the 
engineer. A sample problem including the effect of external stores is 
presented which illustrates in detail the procedures involved in solving 
for the additional root angle of attack, wing load, and other related 
quantities per unit of airplane load factor. 

In order to maintain continuity in the development, the method for 
using the downwash tables from reference 2 is given in an appendix; and 
the pertinent portions of the downwash tables have been included in this 
report with slight modification. 
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wing span, ft 

twist-increment coefficient 
chord, ft 

wing-root chord (at airplane center line), ft 
section normal lift-curve slope 

mean aerodynamic chord 

2 

bending stiffness, lb-in. 

_ * , krth 

downwash factor, w 

r 

2 

torsional stiffness, lb-in. 
acceleration due to gravity 
vortex semi span, ft 

interpolating factor in downwash tables (tables I and II ) 

aerodynamic wing lift outboard of station in question, positive 
for up load, lb 
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L es aerodynamic load acting on external store , lb 

Lfc aerodynamic load acting on horizontal tail, lb 

l running aerodynamic wing load, lb/ft 

M Mach number 

aerodynamic bending moment about axis parallel to airplane 
center line through a given wing station, positive for up 
load outboard of station, in-lb 

My aerodynamic pitching moment about axis perpendicular to air- 

plane center line passing through quarter-chord point at a 
given wing station, negative for up load rearward of axis, 
in-lb 

nigg lift-curve slope of external store in presence of wing per 

degree 

mo two-dimensional section lift-curve slope per degree 

m^ lift-curve slope of rigid wing per degree 

N v number of vortices on both wings 

n normal load factor at airplane center of gravity, positive 

when inertia loads are downward 

q dynamic pressure, lb/sq ft 

S wing area, sq ft 

W gross airplane weight, lb 

w downwash velocity, ft /sec 

X,Y longitudinal and lateral body axes, respectively 

x longitudinal coordinate, positive forward of reference line 

passing through wing quarter-chord at airplane center line 

x C g distance from pitching -moment reference axis to airplane center 

of gravity 

x^ distance from pitching-moment reference axis to aerodynamic 

center of horizontal tail 



k 


y 

a 

Of 

“l 


r 

A 


% 

A c/k 

A 3cA 


lateral coordinate , positive outboard of airplane center line 
angle of attack, deg 

final angle of attack including deformation under airload, deg 
initial angle of attack, deg 

change in section angle of attack due to wing deformation under 
airload, deg 

circulation, sq ft/sec 

difference between two values 

dimensionless vortex spanwise station 

sweep angle of wing quarter-chord line, deg 
sweep angle of wing three-quarter-chord line, deg 


A dominant latent root or eigenvalue 

Matrix notation: 


[] 

LJ 



LI 




U 


-i 



Lil 


square matrix 
row matrix 
column matrix 
diagonal matrix 

rectangular matrix, number of rows not necessarily equal to 
number of columns 

inverse matrix 
error matrix 


unit diagonal matrix 
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||*int|i integrating matrix 

[sj matrix of aerodynamic -influence coefficients per foot; obtained 

from lifting-line theory as derived in reference 1 

[S 2 ] structural elasticity matrix relating change in streamwise 

angle-of-attack distribution to airload distribution, as 
derived in reference 1, ft-deg/lb 

Subscripts and superscripts: 
ac aerodynamic center 

add additional 

adj adjusted 

c chord 

cp center of pressure 

d downwash 

es on external store 

f final, at equilibrium 

i initial 

j exponent in power series 

N number of terms 

n quantity per unit of airplane load factor 

0 opposite wings 

p parabolic 

r at wing root (airplane center line) 

S same wing 

sta station 
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T trapezoidal 

v vortex 

y quantity at semispan station 

a quantity per unit of change in root angle of attack 

Primes indicate the next smaller value. 

METHOD 


The equation relating the normal forces which act on a rigid air- 
plane to produce a change in normal load factor, with wing loads due to 
pitching velocity neglected, is as follows: 

AnW = qm^S Ao. r + AL^ (l) 

When the wing is flexible, that part of the inertia load which acts 
on the wing will cause a change in wing angle— of-attack distribution and 
a resulting increment in aerodynamic lift which is proportional to An. 
Also, the initial aerodynamic load produced by Aoj. will cause the wing 

to deflect and twist with a resulting increment in aerodynamic load which 
is proportional to AOj.. An equation which is equivalent to equation (l) 
but which applies to a flexible wing with negligible pitching accelera- 
tion can, therefore, be written as 

AnW = qm^S Aa^. + AL a Actj. + AI^ An + AL^ (2) 


where at^ is the airload increment due to that part of the wing twist 
caused by a unit change in a*, and ALn is the airload increment caused 
by the wing deflection under an inertia lead of 1 g. 

The quantities AT U and Al^ can be calculated separately. By 
applying the equations of equilibrium between external load and internal 
restraint, the final angle of attack and eerodynamic load produced by 

any initial-angle-of -attack distribution | a i} can be determined. 


Equilibrium-Angle -of -Attach . Distributions 

From reference 1, if the final (at equilibrium) angle of attack "j cu^ 
resulting from the application of |aij- were known, the aerodynamic lift 
distribution could be calculated by the equation 
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For this report, 


{ l f} = [ S l] 1 [ 4c l m o] {“f} 
aA will be defined by 


W = H + M 


(3) 


(4) 


where |a s j is the change in section-angle-of -attack distribution caused 

by the load distribution acting on the flexible wing. As in ref- 

erence 1, f a s} is defined by 


{ a s} - [ s 2] j^f] 


(5) 


Substituting <jaq| + | a s^ f° r e q uat i° n (3) and substituting 

the resulting expression for |z f j, in equation (5) gives 


{ a s} = [ s 2] [ S l] _1 L 4 ^ol ( a i + a s} 


( 6 ) 


Before continuing with the solution for ja s j, it will be desirable 
to factor both the dynamic pressure and compressibility variations out 
of the matrix [^qm^ and convert it to a matrix of constants. This 
operation is performed by multiplying and dividing the matrix by qm R , 
where ir^ is the overall lift-curve slope of the rigid wing. As a 
result. 


[ 4 qmo] = 




mo 

% 


(7) 


mo' 


“R 


can be considered constant at 


Each element of the matrix 
subsonic speeds. Equation (6) can, therefore, be written as 

-li 


[a s ] = qmp [S 2 ] [sj 


m. 


4 — , 

“r! 


| a i + a s} 


( 8 ) 
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Now, solving equation (8) for 



gives 



[I] - qnk[S 2 ] [Sj 


-1 


4 — 


-1 




(9) 


The normal method of solution of equation (9) would be by matrix 


inversion. 


Since the matrix to be inverted 


L 1 ! " <l m R [ S 2] fl] ^ 


is a function of qmp, one inversion must be done for each value of qm^ 

to be investigated. Some method of iteration would be preferable, and a 
very useful method can be developed as follows: By the process usually 

known as "long division," equation (9) can be expanded into a power series 
of any desired number of terms and a remainder term. For convenience, 


let the matrix 
equation (9) may 



be expanded to 


be abbreviated as 



Then, 



(_l] + qn^[S 3 ] + (qn^) 2 ^] 2 + 


+ (qm R ) J [[I] - 



( 10 ) 


where the jth-power term is the exact remainder after J steps in the 
division. In this respect, equation (10) Is analogous to the expansion 


— - — = (l - x) -1 = 1 + x + x 2 + . . . . + (1- x) _1 x J (11) 

1 - X 

Equation (ll) is true for all values of x and J, and equation (10) 
is likewise true for all values of j. Therefore, the answer obtained 
after j steps in the division must be equal to the answer after 
j + 1 steps. In equation form, this may oe stated as 
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Lil + qn R[ S 3] + • • • + (<1 "r) J [W - 1“R [ s 5 j] [s 3 ] J 

J 

= L 1 ] + 1 %[ s 3 ] + • • • + (qm R ) J p J j J + 

(^) J+1 [W - <Pr[S 3 ]] " 1 [S j ] J+1 (12) 

Equation (12) reduces to 

(^R) J [Lll " qm R [S 3 jj [SjJ J 

= ( qm R) J [ S 3] J + ( qm R) J+1 [L I 1 - qm R [ s 3 ]] 1 [ S 3 ] J+1 (15) 

By solving equation (l 3 ) for j" |_f| - qn^js^jj in terms of J and 

P?) ^ ^ shown that, under certain conditions, it will be 

possible to substitute an easily determined scalar quantity for the matrix 

[1/1 - ^RpJ] for the purpose of calculating the remainder term 
(jth -power term) in equation (10). Thus, equation (l 3 ) can be written as 

[W - im R [Sj]J _1 ■= pj] J [S } ]J - qm R p 3 ] J+I j 1 (14) 

When j is sufficiently large, jS^j J+1 is given to a high degree 
of accuracy by Aj^Sjj^, where A is the dominant latent root or eigen- 
value of [Sj] . (See ref. 6 for further discussion.) Substituting 
for [ S 3 J in equation (14) gives 
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[fl - qrDj^Sjj = [sJ J _[S 3 ]J(1 - Aqn^) 


( 15 ) 


However, since 


_ S 5] J ( 1 " ***) 


-1 



l 

1 - Aqm R 


and 



= Id 


equation ( 15 ) reduces to 


W - 1 = W 1 Aqa* 


(16) 


The jth-power term (the remainder) in equation (10) can now be 
written as 


.1 -1 


(q m R ) L 


L 1 ! - ^"r N M 


j .. 


(** r ) 

L - Aqm^ 


w N' 


(i m R) J r= i j 

rr^l'sJ 


Equation (10) in its entirety becomes 



Li] + 'ffv[ s 3] + (^) 2 M 2 + • 


C lm R) 

1 - Aqn^ 


H 


lffi R[ S 3]( a i} 


(IT) 
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It must be made clear that equation (l6) applies only to operations 
on jth-power terms or higher powers of * This does not appear to be 

a practical limitation in aeroelasticity calculations, since increasing 

powers of the typical [^2] [®lj ^ j-~-j matrix rapidly approach the 

L L ^j] 

required simple form. Equation (17) simplifies to 

M - 1%[Sj] {oij + (qmR )2p 3 ] 2 |a 1 j + . . . + [S } ] J+1 |a 1 | 

( 18 ) 

If 

[ S 3]{ a i } = { c e,i} 

and 

M {“i} = { c e,2} = [ s 3]{ c e,4 

then, in general, 

M N { ai } = { c 0,n} = [ s 3]{ c e,N-i} 

It is seen that it is not necessary to raise to the various powers 

if | a ± | is specified. For a specific distribution, equation (18) 

becomes 

W ■ + ( c e,2](^) 2 + ' + { c e,»} 

where N in equation (19) corresponds to J + 1 in equation (l8) . 

In order to determine the number of terms N in equation (19) that 
be required to give any specified accuracy for a trial-and- 

error procedure is used. Essentially, this procedure consists of 
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increasing N until two successive values are found which give the same 
answer to the desired accuracy. In other words, the equation with N 
terms must yield the came result as the equation with N + 1 terms. This 
condition is expressed as 




+ 



1 - Aquip 



+ 


c e,2K<i m R) 2 + 


/ \ N+1 

{ C 0,»}(1"b) N + { C 9,N + l} + { E } (20) 


where 
to /a 



is the difference between the two successive approximations 
Equation (20) reduces to 



and equation (21) can be written as 

{' c e ,N l ' = K4 + 777 W (22) 


Equation (22) represents a system of simultaneous equations with unknowns 

and {e}. The number of equations in the system is equal to the 

1 - Aqmp 1 J 

number of elements in {e|. Therefore, there is one more unknown than 
there are equations. The required extra equation can be obtained in two 
ways: (l) by assuming the value of an element of {e} at a particular 

semispan station to be zero, and (2) by imposing the condition that the 
sum of the squares of the errors must be a minimum (least-squares method). 
If the error is made to be zero at about tte two-thirds-semispan station, 
the first method will give about the same results as the least-squares 

method and will be much simpler. The value of which will give 

1 - Aqmp 

zero error in equation (22) for one semispan station is 
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8 ,N 


'0,N " ^ C 0 ,M 


1 

C 9,K+1 

C 0,N 


qi“R 


(23) 


where C Q N and 

matrices {c 0 ^j 

semispan stations 
Thus, 


< ”'9,N+1 are corres P°nding elements from the column 

and jig ^ ■ ) respectively. The errors at the other 
are easily checked by solving equation ( 22 ) for {e]-. 


r 



) 


The elements of {e j. are then evaluated by inserting the value of 

obtained from equation ( 23 ) into equation (24). It is apparent 


- A<im R 


that 


f. ^ and the error {e} both vary with qn^. Therefore, in 
the trial-and -error procedure the highest value of qn^ to be analyzed 
should be used to calculate - 


Aqn^ 


and the error |e| . 


The error in equation (24) is not the actual error in the 

calculated value of except when {e} s {o}. It can be demonstrated 

that the largest of the errors {e| is greater than the largest true 
error in the assumed solution for ja s j in all cases where the |c q | 

terms alternate in sign. 


When the columns jc^J and + l| have been found which 

ma ^ e sufficiently small, {a g \ can be evaluated from equation ( 19 ) 

which is here repeated: 


| a sJ - j C 0 ,l} qm R + | C 0,2 »(®"r y 


+ 


i 0 > N J1 - AqmR 
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Each initial angle of attack will require a different equation for 
defining the resulting |a s j. As ml 6ht be expected, the number of terms 

required to give any specified accuracy will depend on the shape of the 
initial-angle— of —attack distribution. It can be assumed that the initial- 
angle -of -attack distribution which consists of a unit angle of attack at 

each semispan station will require at least as many terms in its equation 
for ja s | as required by any other distribution normally encountered. 

No harm will be done by carrying more than the minimum number of terms. 


Therefore, if N and — are determined first for the unit sym- 

* ' 1 - Aqmp 

metrical initial-angle-of-attack distribution, it will be unnecessary 
to repeat the trial-and-error process for any other initial angle of 

attack. It will be necessary only to compute the -jCgj. terms through 

/c_ rT l. The same value of — applies to all equations for |<x s !■ 

\ QjNf 1 _ kqmp *• > 


since A is a property of the matrix 



The final-angle-of -attack distribution resulting from a particular 
initial distribution is given by equation (4): 


j^fj- ~ [ a i| + | a s | 


Thus, by combining equations (4) and ( 19) the f inal-angle-of-attack dis- 
tribution is 



+ 



(ggR f_ 

1 - Aqm^ 


(25) 


Airload Distributions 

The airload distribution resulting f rom j a^|- is calculated by 

substituting j. from equation (25) in^o equation (3) and then by 

making use of the previously described conversion of |_4qm 0 "| into 

m n 

amD k _2 to obtain the following equation: 
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(26) 

By performing the term-by-term multiplication indicated in equation ( 26 ), 
an equation for load distribution in terms of qm^ is obtained in the 

form 



and so on for all the terms. 

By integrating equation (27) term by term, an equation for aero- 
dynamic shear (or total lift carried outboard of each station) resulting 

from jap j. is obtained. The integration is done by premultiplying by 

a suitable integrating matrix, which is denoted here by j|l int ||. The 

general form of the shear equation is 



and so on for all the terms. 
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Now, the shear equation is integrated term by term along the Y-axis 
to obtain a bending- (or rolling-) moment equation and along the X-axis 
to obtain a pitching -moment equation. The same integrating matrix which 
was used to obtain shears may be used for the bending moments. If the 
wing quarter-chord locus is straight, the pitching -moment equation can 
be obtained by multiplying the bending -moment equation by -tan A c y 1+ - 

In other cases, a different integrating matrix will have to be devised. 


Application to Calculation of Additional Loading 

The results of the calculations up to this point are equations by 
means of which the final angle of attack, airload, shear, bending moment, 
and pitching moment resulting from a given initial angle of attack can 
be evaluated at any of the chosen wing stations for any value of qm R . 

A particular flight condition may involve several different initial- 
angle -of -attack distributions. Suppose, for example, that it is desired 
to calculate the additional wing loading required to produce a certain 
change in airplane load factor An. The lift forces acting on the air- 
plane are related to the inertia forces by equation (2) 


AnW = qmpS A24. + AI^ Act*. + AI^ An + ALp 


or by 


AnW = ^qm^S + Al^jAOj. + AL n An + AL^ 


The quantity qmp^S + AT ^ is the total wing lift resulting from a 

unit change in wing angle of attack alone j that is, it is twice the aero- 
dynamic shear at the center line calculated from equation (28) when 

is equal to 1 ° at each semispan station. The quantity ALn is 

twice the aerodynamic shear at the center line calculated when |ap}- 

is the twist due to an inertia loading of 1 g on the wing. Equa- 
tion (2) may be written as 


AnW = 2 L r AOj. + 21 ^ An + ALj- ( 29 ) 

where 1^ r and r are the values of Lp at the center line 

calculated by equation (28) for the unit symmetrical and inertial 
twist initial angles of attack, respectively. 
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If a certain value of An is given, the quantities Ao^ . and AL^. 

remain to be determined in equation (29). The additional equation neces- 
sary for this determination is the relationship between the aerodynamic 
and inertial pitching moments acting on the airplane. A convenient ref- 
erence axis for the pitching moments will be that passing through the 
intersection of the wing quarter-chord line with the airplane center 
line. The aerodynamic pitching -moment components corresponding to each 
initial angle of attack will already have been calculated with respect 
to this axis by integrating the shear distributions along the X-axis as 
previously described. The required pitching-moment equation is, thus, 


AnWx cg = 2M Y ^ r AOj, + 34 Y ,n,r An + AL t x t (30) 


where x cg is the distance from the pitching -moment reference axis to 
the airplane center of gravity and x t is the distance from the pitching- 
moment reference axis to the horizontal-tail aerodynamic center. 

By simultaneous solution of equations (29) and (30), Ao^, and at.^ 
are obtained, respectively, as 


An 


W 1 - 


Aop = 


^)^%,n,r--n,r 


2I tc,: 


ST M Y,a,r 


and 


ALt = AnW - 2I^ r - 2L^ r An 


(31) 


(32) 


The additional aerodynamic -shear distribution corresponding 

to the change in load factor An can now be evaluated from its components 
by using the value of A cLp from equation (31). Thus, 


{ L add} = j L a}^ a r + {^An 


(33) 


Si m ilar equations are used to evaluate distributions of angle of attack, 
airload, bending moment, and pitching moment. 
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Antisymmetrical Loadings 


The same principles which have been applied in the previous sections 
to symmetrical loadings may also be applied to antisymmetrical loadings. 

The £s£J matrix for antisymmetrical loadings is different, as will be 

described in the appendix. For this reason, the corresponding matrix 


M N 


4 — 
m R 


is different and a different value of 


obtained for the antisymmetrical case. 


A will be 


Effect of Changing Wing Stiffness 

Investigating the effects of changing the wing stiffness by a cer- 
tain percentage which is constant along the span is a relatively simple 

matter. Such changes have the effect of multiplying [S 2 J by a constant. 

Each term in the angle-of -attack equations needs only to be multiplied 
by the appropriate power of the stiffness ratio. The changes are easily 
followed through the airload, shear, and pitching- and bending -moment 
equations . 


External-Store Effects 


If the wing carries external stores which produce some lift and 


pitching moment, the matrix 




calculated for the wing 


alone is easily adjusted to account for the stores with the following 
provisions : 

(1) The stores should not produce appreciable changes in the addi- 
tional lift distribution of the rigid wing itself. 

(2) The struts supporting the stores should be essentially rigid. 


( 3 ) The variations of store lift and pitching moment with angle of 
attack must be known. 


Each store can be handled separately by the following procedure: 

The lift-curve slope of the external store in the presence of the wing, 
which will be called m„„, can be assumed zo be a constant fraction of 

0 S 

IDp g 

the rigid-wing lift-curve slope m^ so that is constant. The air- 

“R 

load on the store will then be expressed by 
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in 


I S 

'es = a es 


( 34 ) 


from which 


^es %s S 

d ^es ~ qmR ™R~ 2 


(55) 


The airload L es acting at the store center of pressure will produce 

a change in angle of attack at each wing station because of wing twist of 
the flexible wing. The change in wing angle of attack at each station for 
each unit of store load is calculated by the methods described in refer- 
ence 1. The results are arranged in a column matrix, each element of 

da 

which is of the form — — . If each element of the column is multiplied 

dL es . 63 da 

t>y > the resultant column will consist of elements of the form — . 


da. 


es 


da 


Thus, equation (55) can be written as 

da, 


'es 


t-=- = qm R — § es b 

da es K d 


“es 


“R 2 


(56) 


or 


1 da s 


qm^ da^ 


da s m es S 
-es ^es “R 2 

The typical element of the column matrix 


has the same prop- 


1 [ da s 

^Ri^esJ 

erties as an element of the matrix j~J^S . Under the assunp- 

tion that the external-store strut is rigid, a„„ will be the same as 

6S 

the wing angle of attack at the store semispan station. If the store 
station coincides with one of the vortex stations, the column matrix 

can be added directly to the column of the matrix 


1 

da s 

qn^ 

da es 
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M [ S J 


-1 


4 ^ 

m R 


which corresponds to that vortex station. In the 


general case, where the store is between two stations, the column matrix 
da 0 


qm R 


da. 


is divided between the two corresponding columns so that 


es 


the greater portion is given to the column of 


ft * 1 

L 


nr 


4 

“r 


corresponding to the closer vortex station. For instance, if the store 
is located at a point one-third of the distance between two stations, 


two-thirds of the value of 


_iK 

<Pr { da esJ 


is added to the nearest adjacent 


column and the remainder is added to the other adjacent column of 




This procedure is equivalent to Interpolating 


between two wing angles of attack to determine the store angle of attack 
When more than one store is present, several columns of 

[ S 2] [ S l] “ 1 [ 4 wil1 be affected. 

The jCgj columns and all quantities derived from them must be 
recalculated by using the adjusted 


NN't 2*1 


matrix. After 


the equation for has heen determined, an interpolation is made 

between the column elements corresponding to the two wing stations adja- 
cent to each store to allow evaluation of the final angle of attack at 
the store stations for any value of qm^. The load acting on each exter 

nal store at equilibrium is then 


^es 


m 


es S 


^ ,es qm R i“~ 2 


(37) 


This load must be added to the new shear equations at each station 
inboard of the store. Each station inboard of the store is also sub- 
jected to a bending moment and pitching moment caused by the load acting 
on each store. These moment increments are 


^b,es “ L es ^es 
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and 


^Y,es ^es ^es 

where Ay es and Axgg are measured from the external- store center of 
pressure to the wing quarter-chord point at each vortex station involved. 

APPLICATION OF METHOD TO SAMPLE PROBLEM 


Sample calculations of the additional wing loading have been worked 
out for the wing-nacelle combination shown in figure 1. The calculations 
are divided into three parts. The first part, entitled "influence Coef- 
ficients," includes the determination of the , ^SgJ , and |^4 — 

matrices and of the external- store effects. The second part, entitled 
Load-Distribution Coefficients," includes the determination of equa- 
tions for the required components of final angle of attack, airload, 
shear, bending moment, and pitching moment in terms of qmp . The third 

P ai> t, entitled Additional Airloads," consists of the determination of 
the wing-root airloads and tail airloads per unit of load factor for 
specific conditions of qm R , gross weight, and center-of -gravity location. 

For the sample problem, the lift of the wing was represented by 
20 horseshoe vortices, each 10 percent of the semispan in width. (See 
fig. 1.) The center of each of the bound vortices is located on the 
wing quarter-chord locus. The El and GJ distributions which were 

used to calculate the pSgj matrix by the method of reference 1 are 

plotted against nondimens ional semispan station in figure 2. The dimen- 
sions of the sample wing are as follows: 


Span, ft ± 16 

Vortex semi span, h, ft 2.9 

Area, sq ft 1,428 

Root chord, ft 17-34 

Taper ratio 0.42 

Aspect ratio 9.42 

Quarter -chord sweep angle, deg 55 

Three-quarter-chord sweep angle, deg 31-5 

External- store location, percent of semispan 38.2 
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Influence Coefficients 


The inverse of the aerodynamic-influence-coefficient matrix 


the matrix 4 — , and the structural-influence-coefficient matrix 

L “r 

js 2 "{ are required in the calculation of final airload distribution, as 
shown by equation (26): 

•M = ^N" 1 4 ^ 'W + qm R{ C 9,l} + • • • + *qmJ C0 > N }{ 

v ^ 


where 




| C 0,n} = [ S 3]{ C 0,N-l} 

The product [S 2 ] jsj -1 |4 ^jj is first calculated for the wing alone, 

and then the method of adding external- store effects is demonstrated. 

Calculation of [Sj .- The method of deriving the £s J matrix 
from the downwash tables (tables I and II, is described in detail in 
the appendix. For the sample problem, the required j^j matrices 

are those given in equations (a 6) and (A7, ; • The necessary matrix 

is calculated from equation (All) as 


tan 31. 5° + IjX; - tan 35° 


which becomes, numerically, 
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fir]- 


1.3421 2.7425 

4.1429 5.5434 

6.9438 

8.3442 

9 . 7446 

11.1450 

12.5455 

13.9459 

0.1150 1.5154 

2.9158 4.3162 

5.7167 

7*1171 

8.5175 

9.9179 

11.3183 

12.7188 

-1.1121 0.2883 

1.6887 3.O891 

4.4895 

5.8900 

7.2904 

8.6908 

10.0912 

11.4916 

-2.2763 -0.9588 

0.4616 1.8620 

3*2624 

4.6628 

6.0633 

7*4637 

8.8641 

10.2645 

-3.5664 -2.1660 

-0.7655 0.6349 

2.0353 

3*4357 

4.8361 

6.2366 

7.6370 

9.0374 

-4.7935 -3.3931 

-I.9927 -0.5922 

0.8082 

2.2086 

3.6090 

5.0094 

6.4099 

7*8103 

-6.0206 -4.6202 

-3.2198 -1.8194 

-0.4189 

O.9815 

2.3819 

3.7823 

5.1827 

6.5832 

-7.2477 -5*8473 

-4.4469 -3*0465 

-1.6461 

-0.2456 

1.1548 

2*5552 

3.9556 

5.3560 

-8A749 -7.0744 

-5.6740 -4.2736 

-2.8732 

-1.4728 

-0.0723 

1.3281 

2.7285 

4.1289 

-9*7020 -8.3016 

-6.9011 -5*5007 

-4.1003 

-2.6999 

-1.2995 

0.1010 

1.5014 

2.9018 


(38) 


By the use of 


*3L 

h 


from equation (38) 


from equation (a 6) and 

as arguments in the downwash tables (tables I and II ), the downwash fac- 
tors for vortex and downwash stations on the same wing ^FgJ are deter- 


mined to be as follows: 


['■]- 


4.4942 

-1.2369 

-0.2313 

-O.O963 

-0.0527 -0.0332 

-0.0228 -O.OI63 

-0.0127 

-0.0100 

-0.7176 

4.3964 

-1.2454 

-0.2331 

-O.O969 -0.0530 

-0.0334 -0.0229 

-0.0163 

-0.0128 

-0.0955 

-0.7919 

4.3244 

-1*2530 

-0 . 2349 -0 . 0976 

-0.0532 -0.0335 

-0.0230 

-0.0164 

-0.0364 

-0.1010 

-0.8608 

4.2835 

-1.2596 -0.2364 

-0.0982 -0.0535 

-0.0337 

-0.0231 

-O.OI87 

-0.0373*- 

-O.IO68 

-0.9224 

4.2286 -1.2656 

-0.2382 -O.O987 

-0.0538 

-0.0338 

-0 . 0114 

-0.0192 

-0.0387 

-0.1125 

-O.9761 4.1955 

-1.2712 -0.2393 

-0.0992 

-0.0540 

-0.0077 

-0.0117 

-O.OI98 

-0.0401 

-O.H7I -1.0222 

4.1692 -I.2756 

-0.2404 

-0 . 0997 

-0.0057 

-0.0078 

-0.0120 

-0.0204 

-0 . 04l6 -0.1246 

-1.06l4 4.1479 

-1.2798 

-0.2415 

-0.0042 

-0.0057 

-0.0080 

-0.0122 

-0.0209 -0.0431 

-0.1307 -1.0949 

4.1301 

-1*2833 

-0.0032 

-0.0042 

-0.0057 

-0.0081 

-0.0125 -0.0215 

-0.0447 -0.1368 

-1.1234 

4.1154 


Similarly, by using 


Ay 
h J 


in place of 


[f| . 

L h J 0 


(39) 


the downwash factors 


L “ J 0 LI1J S 

for vortex and downwash stations on opposite wings are found to be 
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-0.0014 

-0.0017 

-0.0019 

-0.0025 -0.0027 

-0.0035 

-0.0040 

-0.0049 

-0.0064 

-0.0016 

-0.0018 

-0.0021 

-0.0025 -0.0051 

-0.0057 

-0.0047 

-0.0059 

-0.0075 

-0.0017 

-0.0020 

-0.0024 

-0.0026 -0.0055 

-0.0045 

-0.0054 

-0.0071 

-0.0092 

-0.0018 

-0.0022 

-0.0026 

-0.0052 -0.0040 

-0.0050 

-0.0065 

-0.0085 

-0.0117 

-0.0020 

-0.0024 

-0.0029 

-0.0056 -0.0045 

-0.0060 

-0.0079 

-0.0107 

-0.0152 

-0.0021 

-0.0026 

-0.0052 

-o.oo4o -0.0052 

-0.0069 

-0.0095 

-0.0157 

-0.0206 

-0.0025 

-0.0028 

-0.0056 

-0.0046 -0.0061 

-0.0085 

-0.0121 

-0.0181 

-0.0295 

-0.0025 

-0.0051 

-0.0059 

-0.0052 -0.0071 

-0.0101 

-0.0154 

-0.0252 

-0.0461 

-0.0026 

-0.0054 

-0.0045 

-0.0058 -0.0082 

-0.0125 

-0.0200 

-0.0569 

-0.0814 

-0.0026 

-0.0056 

-0.0047 

-0.0065 -0.0094 

-0.0152 

-0.0266 

-0.0581 

-0.1826 


-0.0079 

-0.0096 

-0.0125 

-0.0165 

-0.0224 

-0.0527 

-0.0521 

-0.0956 

-0.25H 

-1.2448 


The aerodynamic- influence -coefficient matrix for symmetrical 

airloads, which is found by dividing the vortex semispan h into the 
sum of the downwash -factor matrices (eqs. (39) and (^0)), is written as 

N - [N ♦ N] i 

or, numerically, as 


1.5492 
- 0 . 2480 
-0.0555 
-0.0152 
-0.0071 
-0.0047 
-0.0054 
-0.0028 
-0.0025 
- 0.0021 


-0.4271 

-0.0804 

-0.0540 

-0.0191 

-0.0126 

-0.0092 

-0.0073 

-0.0066 

1.513* 

-0.4502 

-0.0812 

^0.0545 

-0.0196 

-0.0151 

-0.0099 

-0.0082 

-0.2758 

1.4904 

-0.4550 

-0.0822 

-0.0351 

-0.0202 

-0.0140 

-0.0111 

-0.0556 

-0.2977 

1.4760 

-0.4357 

-0.0^52 

-0.0561 

-0.0214 

-0.0157 

-0.0157 

-0.0578 

-0.5193 

1.4566 

-0.4J95 

-0.0849 

-0.0377 

-0.0258 

-0.0075 

-0.0144 

-0.0402 

-0.3584 

1.4144 

-0.4416 

-0.0872 

-0.0415 

-0.0050 

-0.0081 

-0.0154 

-0.0425 

-0.5353 

1.4555 

-0.4461 

-0.0951 

-0.0058 

-0.0055 

-0.0088 

-0.0168 

-0.0^64 

-0.5713 

1.4216 

-0.4572 

-0.0051 

-0.0042 

-0.0062 

-0.0100 

-0.03 91 

-0.0520 

-0.3903 

1.5961 

-0.0027 

-0.0056 

-0.0050 

-0.0076 

-0.0: 27 

-0.0246 

-0.0672 

-0.4505 

on of 

4 ?a| 

m R 

The 

N 

matrix can 

now be used 


-0.0062 

- 0.0078 

- 0.0100 

-0.0156 

-0.0194 

-0.0299 

-0.0525 

-0.1162 

-0.5222 

0.9899 


late the matrix 4 — from wind-tunnel measurements of pressure dis- 

®R 

tribution on rigid-wing models. This process is described in detail in 
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reference 1, but, briefly, from the definition of j^S-[| , the following 
equation is derived: 


[ s i]{ i }‘ MM 


When the model is rigid, any change in angle of attack can be 
assumed constant along the semispan. Taking the derivative with respect 

to a and factoring q out of |4qmo| , thus, results in 
or 

MW-W 

However, 


dll 1 
da < 1 



the elements of which are easily calculated from the chordwise pressure 
distributions at the required semispan stations. ^The chords c must 

be those of the full-scale airplane to agree with [s£] .) Therefore, 


[ s i]{v} = Nw - H} 


(45) 


Dividing equation (43) by the wing lift-curve slope m R gives 


1_ 

m R 




The lift-curve slope 


“r 


is calculated from 



“R = 


Integrating row 


j{° n a C } 


2 

S 


(44) 


(45) 
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where the row matrix is taken frcm the bottom row of the integrating 
matrix j| ^ int >p {{ 6^ ven table HI and is multiplied by the wing semi- 
span given in feet. 


The measured values of 



for the sample wing are as follows: 



0.4746 
0.7034 
0.8256 
0.9155 
+ 0.9829 „ 
1.0402 
1.0943 
1.1364 
1.1702 
1.2001 


The value of m^ from equation (45) is, thus, O.O768I. 


The diagonal matrix 


as 


4 

m R 


is then calculated from equation (44) 


5.5391 0 

0 5.4665 

0 0 

0 0 



0 0 
0 0 
0 0 
0 0 


0 

0 

5.7976 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

5.9967 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

5.8882 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

5.8x59 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

5.9028 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

5.8988 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

6.0459 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

6.8564 


M 


( 46 ) 
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Calculation of MS 


In the calculation of 


-1 4 ^ 


9 the aerodynamic -inf luence -coefficient matrix 


inverted and postmultiplied by 


4 ^ 

_ “r 


with the result that 




2.4208 1.1850 0.6784 0.4565 0.5274 0.2560 0.2208 0.2005 0.2002 0.2554 

0.4507 4.0852 1.4780 0.8570 O.5662 0.4195 0.5478 0.5070 0.5007 0.5478 

0.1671 0.8975 4.5407 1.6797 0.9592 0.6555 0.5148 0.4579 0.4190 0.4786 

0.0860 0.3586 1.1221 4.8695 1.7919 1.0624 0.7717 0.6240 0.5791 0.6520 

0.0543 0.1929 0.4782 1.3535 4.9754 I.9229 1.2178 0.9142 0.8120 0.8964 

0.0392 0.1251 0.2660 0.5998 1.4858 5.0978 2.1242 1.3974 1.1614 1.2457 

0.0312 0.0924 0.1772 0.3468 O.7073 1.6532 5.3869 2.3547 1.7317 1.7722 

0.0266 0.0750 0.1343 0.2394 0.4312 0.8460 1.9104 5-6812 2.8230 2.6282 

| 0.0235 0.0650 0.1115 O.I878 0.3134 0.5542 1.0766 2.2707 6.3971 4.2354 

^0.0216 0.0587 O.O989 0.1623 0.2611 0.4400 0.7970 1.5081 3.1719 9.1053 


Calculation of 




-1 , m Jl 

4 — .- In the calculation of 




} the structural elasticity matrix 


by using the method of reference 1 is 


-1 . 2282 

-1.1817 

-1.0408 

-0.8509 

-0.6538 

-0.4739 

-0.3234 

-0.1979 

-0.08814 

0 

-0.7709 

-0.7772 

-0.7503 

-O.6561 

-O.5258 

-0.3923 

-0.2732 

-0.1696 

-O.076H 

0 

-0.4739 

-0.4739 

-0.4786 

-0.4613 

-0.3978 

-0.3107 

-0.2232 

-0.1413 

-0.06404 

0 

-0.2763 

-O.2763 

-0.2763 

-0.2801 

-O.2699 

-0.2292 

-0.1731 

-0.1130 

-0.05201 

0 

-0.1488 

-0.1488 

-0.1488 

-0.1488 

-0.1521 

-0.1476 

-0.1230 

-0.08471 

-0.03998 

0 

-0.07082 

-0.07082 

-0.07082 

-0.07082 

-0.07082 -O.O7395 

-0.07288 

-0.05640 

-0.02792 

0 

-0.02596 

-0.02596 

-0.02596 

-O.O2596 

-0.02596 -0.02596 

-0.02858 

-0.02808 

-0.01589 

0 

-O.OO236 

-O.OO236 

-0.00236 

-0.00236 

-0.00236 -0.00236 

-0.00236 

-0.00439 

-0.00386 

0 

0.00302 

0.00302 

0.00302 

0.00302 

0.00302 

0.00302 

0.00302 

0.00302 

0.00150 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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The product of equation (47) premultiplied by equation (48) gives 


-0.3435 

-0.5167 

-0.4529 

-0.3612 

-0.2652 

-0.1932 

-0.1469 

-0.1166 

-0.1072 -0.1222 

-0.3526 

-0.5158 

-0.4507 

-0.3597 

-0.2640 

-C.I923 

-0.1461 

-0.1158 

-0.1064 -0.1214 

-0.2975 

-0.4865 

-0.4392 

-0.3517 

-0.2583 

-0.1878 

-0.1423 

-0.1123 

-0.1030 -0.1174 

-0.2468 

-0.4241 

-0.4050 

-0.3339 

-0.2458 

-0.1783 

-0.1343 

-0.1051 

-0.0958 -0.1091 

-0.1921 

-0 . 3416 

-0.5426 

-0.2986 

-0.2258 

-C.1632 

-0.1217 

-0.0941 

-0.0850 -0.0966 

-0.1407 

-0.2564 

-0.2664 

-0.2441 

-0.1949 

-0.1437 

-0.1056 

-0.0800 

-0.0714 -0.0809 

-0.0969 

-0.1797 

-0.1914 

-0.1817 

-0.1526 

-0.1183 

-0.0875 

-0.0642 

-0.0562 -0.0635 

-0.0597 

-0.1121 

-0.1214 

-0.1180 

-0.1025 

-C.O834 

-0.0642 

-0.0464 

-0.0391 -0.0441 

-0.0267 

-0.0504 

-0.0551 

-0.0543 

-0.0480 

-0.0400 

-0.0318 

-0.0233 

-0.0191 -0.0214 


External- store effects .- The presence of an external store requires 

the modification of the matrix equation = |[s 2 ] js£J 1 |4 ^|j as 

described in the "Method" section. The external store at the 58.2-percent- 
semispan station is estimated, from wind-tunnel data, to carry about 

2 percent of the rigid-wing airload ^ - = 0 . 02 ^ from which the wing 

twist due to store load per unit of store angle of attack is calculated 
as 

0.0005851 
0.0005851 
. 0.0005851 

da s _ J 0.00058511 (50) 

Idoesf |0- 0005851 

< J 0.0005651 

0.0004710 
0.0002381 
0.0000791 
l 0 

The store is located between vortex stations 6 and 7 = 0.45 

and O. 55 ). The store angle of attack, which is assumed to be the same 
as that of the wing at the store station, is by linear interpolation 


0^3 = 0 . 52 Og + 0.680^ 


If 32 percent of 


is added to column 6 of the 
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(eq. (4-9)) and. the remainder is added to column 7> "the external- store 
effect on wing twist will be included when the adjusted matrix ^ 


which is given as 


-0.3435 -0.5167 

-0.4529 

-0.3612 

-0.2652 

-0.1745 

-0.1071 

-0.1166 

-0.1072 

-0.1222 

-0.3326 -0.5138 

-0.4507 

-0.3597 

-0.2640 

-0.1736 

-0.1063 

-0.1158 

-0.1064 

-0.1214 

-0.2973 -0.4865 

-0.4392 

-0.3517 

-0.2583 

-O.I69I 

-0.1025 

-0.1123 

-0.1030 

-0.1174 

-0.2468 -0.4241 

-0.4050 

-0.3339 

-0.2458 

-0.1596 

-0.0945 

-0.1051 

-0.0958 

-0.1091 

-0.1921 -0.3416 

-0.3426 

-0.2986 

-0.2258 

-0.1445 

-0.0819 

-0.0941 

-0.0850 

-0.0966 

-0.1407 -0.2564 

-0.2664 

-0.2441 

-0.1949 

-0.1250 

-0.0658 

-0.0800 

-0.0714 

-0.0809 

-0.0969 -0.1797 

-0.1914 

-0.1817 

-0.1526 

-0.1032 

-0.0555 

-0.0642 

-0.0562 

-0.0635 

-0.0597 -0.1121 

-0.1214 

-0.1180 

-0.1025 

-0.0764 

-0.0493 

-0.0464 

-0.0391 

-0.0441 

-0.0267 -0.0504 

-0.0551 

-0.0543 

-0.0480 

-0.0374 

-0.0265 

-0.0233 

-0.0191 

-0.0214 

0 0 

0 

0 

0 

0 

0 

0 

0 

0 


is postmultiplied by the column matrix of wing angle of attack (eq. (50)). 


Load-Distribution Coefficients 

The next series of calculations evaluates the coefficients required 
to determine the distribution of final angle of attack, airload, shear, 
bending moment, and pitching moment resulting from (a) a unit change in 
wing angle of attack and (b) the initial wing twist due to wing inertia 
loading per unit of change in load factor. Both final-angle-of -attack 
distributions are calculated by the use of equation (25), which requires 
an iteration procedure to determine the number of terms necessary. 

Equations for final-angle-of-attack distribution .- The coefficients 
in the equation for the final-angle-of-attack distribution (eq. (25)) 
resulting from a unit symmetrical initial -angle-of -attack distribution 
are evaluated in table IV. The column in table TV is, in this case, 

unity at each station. The values under |Cq -A are by definition 


and the succeeding columns are defined as 

{ c e,t>} = [ s 3] ad j{ c e,i«-i} 
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Each { c e} column In table IV, therefore. Is obtained by premulti- 


plying the previous column by S 


adj 


from equation (51) • After three 


w columns were calculated, a check was made to see whether the errors 


resulting from using only the first two 


{ c e> 


columns would be excessive. 


In order to perform the check, a value of was estimated from 

1 - Aqm R 

the values of |Cq 2 ^ an d appearing at the 65 -percent-semispan 

station (station 4-) with qm R taken to be 100 lb y^t^-deg and with the 

use of equation ( 23 ). The value of -■ thus obtained was recorded 

1 - ?\qm R 

at the foot of the |Cq^ 2 | column and was used to obtain the errors in 
the Hj -2 column by applying equation (24). These errors were con- 
sidered too large. Therefore, another column ^j. was computed and 

the check was repeated. This time the errors, listed in the 

column, were considered satisfactorily small. A sufficiently accurate and 
convenient form for the angle-of -attack equations was, therefore, con- 
sidered to be 


{ a f,a} = ( ai j + { C e,lj q?n R + | c e,2}(«“R) 2 + { c 0,3j( qm R ) 5 
Equation (52) may be rewritten with the numerical values of the 


coefficients from table IV as follows: 


1.0 

-0.025671 

0.051241 x 10 -2 

1.0 

-0.025442 

0.050743 X 10 -2 

1.0 

-0.024374 

0.048428 x 10 -2 

1.0 

-0.022198 

0.043760 x 10" 2 

1.0 

-0.019028 

0.037071 x 10~ 2 

1.0 

-0.015256 

0.029291 x 10 -2 

1.0 

-0.011447 

0.021565 x 10 -2 

1.0 

-0.007690 

0.014141 x 10 -2 

1.0 

-0.003622 

0.006550 x 10 -2 

1.0 

0 

0 



1 

1 - Aqmp 

(52) 

{°e} 


- 0 . 101 X 64 x 10" 4 
-0.100176 x 10" 4 
- 0.095584 x 10 -4 
- 0.086353 x 10 -4 

_k 

-0.073092 x 10 
-0.057710 x 10 _4 
- 0.042449 x 10 -4 
-0.027801 x 10 -4 
-0.012868 x 10 -4 
0 


1 

( qm H ) 2 

(^“r ) 3 


1 - 


( 55 ) 
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Equation (53) gives the final angle of attack ja f resulting 

frcm a symmetrical angle -of -attack change of 1°. An equation of similar 
form which gives the final-angle-of -attack distribution <ja f n | caused 

by the wing twisting under an inertia loading of 1 g is as follows: 



0.5208 

-0.0098555 

0.019598 x 10~ 2 

-0.058266 x 10 -4 

0.5094 

-0.0097554 

0.019209 x 10 -2 

-0.057892 x 10" 4 

0.4859 

-0.0095007 

0.018527 x 10' 2 

-0.056155 x 10' 4 

0.4-515 

-O.OO8587I 

0.016541 x 10' 2 

-0.052654 x 10 -4 

0.5558 

-0.0070858 

0.014011 x 10' 2 

-0.027645 x 10" 4 

0.2246 

-0.0055770 

0.011060 x 10" 2 

-0.021826 x 10 _4 

0.1507 

-0.0040856 

0.008154 x 10 -2 

-0.016055 x 10 -4 

0.1529 

-0.0026597 

0.005525 x 10 -2 

-0.010515 x 10 -4 

O.O695 

-0.0012269 

0.002464 x 10 -2 

-0.004866 x 10 -4 

0 

0 

0 

0 


l 

< (« d r) 2 

1 - Aqnip 




(54) 


The first column of the matrix of coefficients in equation (54) is the wing 
angle -of -attack distribution due to dead— weight loading* Each succeeding 

column is the product of the previous column premultiplied by fsjl 

L ^ J adj 

It 1® not necessary to repeat the iteration procedure, because the number 
of coefficient columns required to obtain acceptable accuracy is only 
slightly dependent on the shape of the initial-angle-of -attack 
distribution. 


External-store angle of attack .- In order to calculate the airload 
carried on the external store, the angle of attack at the store loca- 
tion must be known. The final angle of attack of the store resulting 
from a unit change in wing-root angle of attack is obtained by interpo- 
lation between the rows of equation (53) which correspond to stations 6 
and 7. Thus, 


a, 


es,a 



-2 


(Pp 

C 5 ^) 2 > 

(n ) 3 


1 ~ 

( 55 ) 


-0.012666 0.024037 x 10 


-0.047333 x lcf^j 
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Similarly, the following equation for store angle of attack resulting 
fran the unit- load-factor inertia loading is obtained by interpolation 
in the rows of equation (54): 

1 

qmR 


a 


'es 


n = [0.1743 -0.004563 0.009070 x 10" 2 -0.01790 x 10- 4 , (<P R ) 

(«“r) 


1 - Aqm^ 

(56) 


Running-load-distribution equations . - The equation for running-load 
distribution jl f a j» resulting from a unit change in wing-root angle of 
attack is obtained by premultiplying the coefficients of equation (55) 

-1 


by qn^ [sj 


|4 ^2. 

®R 


Thus, 



‘f. 


6.17868 

-0.13389 

O.26554 x 10" 2 

-0.52407 x 10” 4 

9 . 15804 

-0.19511 

0-38646 x 10” 2 

-O.76266 x 10” 4 

10 . 74948 

-0.21411 

0.42228 x 10" 2 

-0.83316 x 10” 4 

11.91720 

-0.21281 

0.41795 x 10” 2 

-0.82211 x 10” 4 

12 . 797*+0 

-0.19548 

0.37936 x 10” 2 

-0.74784 x 10” 4 

13.54236 

-0.16930 

0.32506 x 10 ” 2 

-0.64046 x 10“ 4 

14.25372 

-0.13989 

0.26550 x 10” 2 

-0.52281 x 10 ” 4 

14.795^0 

-0.10807 

0.20308 x 10“ £ 

-0.39942 x 10” 4 

15.23520 

-0.07651 

0.14320 x 10” £ 

-0.28179 x 10” 4 

15 . 62484 

-0.05195 

0.09783 x 10" 2 

-0.19257 x 10” 4 


qmp 

( qm R) 2 

Y 


1 - 


(57) 


i'f.n} 


resulting from the 


The equation for running-load distribution 
twist due to an inertia loading of 1 g is obtained by premultiplying the 

coefficients of equation (54) by qn^pS^ 4 —I • Thus, 
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1.39791 

-0.05098 

0.10048 

X 

icr 2 

1 

0 

H 

£ 

X 10 -4 

3.80110 

-0.07416 

0.14621 

X 

10’ 2 

-0.28848 

X 10" 4 

4.10930 

-0.08095 

0.15920 X 

10‘ 2 

-0.31514 

x 10" 4 

3-95352 

-0.07984 

0.15755 x 

10~ 2 

-0.31094 x 10' 4 

3.39820 

-0.07242 

0.14333 

X 

10’ 2 

-0.28284 

X 10~ 4 

2 . 70251 

-0.06188 

0.12274 

X 

10’ 2 

-0.24222 

x 10 -4 

2.15264 

-0.05039 

0.10018 

X 

10’ 2 

-0.19772 

x 10 -4 

1 . 79402 

-0.03840 

O 

b 

CT\ 

X 

10’ 2 

-0.15105 

X 10 -4 

1.31759 

-0.02707 

0.05398 

X 

10’ 2 

-0.10666 

x 10' 4 

0.88601 

-0.01853 

ON 

CO 

CS 

0 

6 

X 

10’ 2 

-0.07282 

x 10“ 4 


1 - 


Shear-distribution equations.- The matrix of coefficients for the 
shear-distribution equations is calculated by premultiplying the matrix 


of running-load coefficients by the integrating matrix j|lint,p||j w ^ich 

is given in table III. This integrating matrix is based on Simpson's 
parabolic method and is given in nondimensional form for 10 stations, 
with an extra row added to give the integral to the center line. The 
results are made dimensionally correct by multiplying them by the wing 
semispan (in feet, for the sample problem). 

The complete shear equation must also include the concentrated load 
acting on the nacelle. This load which is obtained by substituting 
a es a ^ rom equation (55) into equation (3*0 is given by 


2 


-0.012666 0.024037 x 10" 


-0.047333 x 10-J ^ ^3 

(^r) 4 

1 - Xqm R 


Since = 0.02 and — = 71^ square feet, equation (59) can be written 

m R ^ 

as 

, , ( qm R) 2 


Le S}a = 14.3 -0.18 0.34 x 10~ 2 -0.68 x lO"^ 




3 L (60) 


(^"R) 

1 - Mm R 
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Equation ( 60 ) gives the external-store lift per unit of wing-root angle 
of attack. The store airload due to inertia twist per unit of load fac- 
tor is calculated similarly by substituting ote S>n from equation (56) 

into equation ( 34 ) to obtain 


qmR 

, (®“r ) 2 1 

Les, n = L 2 ' 5 "°'° 7 °' 15 X 10 " 2 “°- 26 X 10 J I (^“ r ) 5 [ ^ 6l) 

(q^R ) 4 


1 - 


Ml 


Since the external-store lift acts as a concentrated load on the 
wing, it is to be added to the shear at each station inboard of the store. 
For the unit symmetrical angle -of -attack distribution, integration of 
equation ( 57 ) by using table III and the addition of the external-store 
coefficients from equation (60) ^with the qm^ terms factored out^ 

result in the following equation for shear distribution: 


N 



Load 

on wing: 


9.7 

-0.21 

0.0042 

-0.000083 

54.9 

-1.19 

0.0235 

-0.000464 

112.8 

-2.38 

0.0471 

-0.000931 

178.7 

-3.65 

0.0717 

-0.001414 

250.4 

-4.82 

0.0949 

-0.001870 

326.8 

-5.88 

0.1153 

-0.002274 

407.5 

-6-77 

0.1325 

-0.002611 

491.8 

-7.49 

0.1460 

-0.002878 

578.9 

-8.02 

0.1560 

-0.003075 

668.4 

-8.39 

0.1629 

-0.003211 

i 714.0 

-8.53 

0.1655 

-0.003262 


Load on external store: 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

14.3 

-0.18 

0.0034 

-0.000068 

14.3 

-0.18 

0.0034 

-0.000068 

1^.3 

-0.18 

0.0034 

-0.000068 

14.3 

-0.18 

0.0034 

-0.000068 

14.3 

-0.18 

0.0034 

-0.000068 


qn^ 
(<J»r ) 2 
' (i“r ) 5 - 
( < P'r ) 4 

1 - Aqm^ 


(62) 


Similarly, the equation for shear distribution due to wing-inertia 
twist per unit of load factor is obtained from equations (58) and (6l) 


as 
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M = 


Load 

on wing: 

2.1 -0.08 

0.0016 -0.000051 

18.1 -0.45 

O.OO89 -0.000175 

41.5 -0.90 

0.0178 -0.000552 

64.9 -1.57 

0.0271 -0.000555 

86.5 -1.82 

0.0558 -0.000707 

105.9 -2.21 

0.0456 -0.000860 

117.9 -2.53 

0.0500 -0.000988 

129.4 -2.79 

0.0551 -0.001089 

158.4 -2.98 

0.0589 -O.OOII65 

144.8 -5.11 

0.0615 -0.001214 

147.0 -5.16 

0.0625 -0.001254 


Load on external store: 


0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


2.5 -0.07 0.0015 -0.000026 

2.5 -0.07 0.0015 -0.000026 
2.5 -0.07 0.0015 -0.000026 

2.5 -0.07 0.0015 -0.000026 

2.5 -0.07 0.0015 -0.000026 




(®“Rr 

ia£ 


Jl - Xqmpj 


(65) 


The shear coefficients due to distributed wing lift and concentrated 
store lift are to be kept as two separate matrices for ease in obtaining 
the bending-moment and pitching-moment equations. 

Bending-moment-distribution equations . - The matrix of coefficients 
for the component of the shear equation due to distributed wing lift is 
Integrated to obtain the matrix of coefficients for the component of 
bending moment due to distributed wing lift. For this purpose the inte- 
grating matrix ||l int ^ p || in table III multiplied by b/2 may be used. 

However, shear-distribution curves are usually straight enough that the 
much simpler trapezoidal type of integrating matrix can be used. This 

integrating matrix ||li n t,T|j is given in table V and is nondimens ion- 

alized with respect to semispan. For use in obtaining bending moments 
in inch-pounds, the elements of table V are multiplied by the semi span 
in inches. 

The coefficients for calculating the bending moments due to the 
concentrated load imposed by the external store are obtained by pre- 
multiplying each row of the matrix of shear coefficients due to concen- 
trated store load by a diagonal matrix whose elements are 
y es - y s t a = Ay es a For the sample problem, 
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By using equation (62), table V, and [Ay es j from equation (64) in 

the prescribed manner, the complete equation for bending moment per unit 
of wing-root angle of attack is given as follows: 


Load on wing: 
338 -7-3 0.145 

2,585 -55-9 1-109 

8,420 -180.0 3.567 

18,563 -389.2 7-702 

33,495 -683.0 13.498 

|Mb,a}- - 53,582 -1,055-0 20.814 

79,136 -1,495.2 29.437 

110,432 -1,991.7 39.130 

147,692 -2,531.6 49.642 
191,098 -3,103-0 60.742 
214,359 -3,595.0 66.412 


Loa 1 on external store: 

-0.00287 000 0 

-0.02188 0000 

-0.07040 000 0 

-0.15198 000 0 

-0.26627 0000 

-0.41048 +000 0 

-0.58047 318 -4.0 0.076 -0.00151 

-0.77151 1,312 -16.6 0.315 -0.00621 

-O.97867 2,306 -29.2 0.554 -0.01091 

-1.19740 3,300 -4l.8 0.793 -0.01562 

-I.30913 3,796 -48.1 0.913 -0.01798 



The equation for bending moment due to inertia twist per unit of 
load factor, obtained from equation ( 63 ) in a fashion similar to that 
of equation ( 65 ), is as follows: 
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Load on wing: 


Load on external store: 


K4 ■ 


5 72 

-2.8 

0.055 

-0.00109 


0 

0 

0 

0 

1 775 

-21.3 

0.420 

-0.00828 


0 

0 

0 

0 

j 2,844 

-68.5 

1.349 

-0.02663 


0 

0 

0 

0 

6,540 

-147.8 

2.912 

-0.05749 


0 

0 

0 

0 

11,800 

-258.8 

5.100 

-0.10071 


0 

0 

0 

0 

18,418 

-398.9 

7.862 

-0.15526 

+ 

0 

0 

0 

0 

26,136 

-563.8 

11.119 

-0.21956 


55 

-1.5 

0.029 

-0.00057 

34,741 

-749-0 

14.778 

-0.29181 


229 

-6.0 

0.119 

-0.00235 

44,060 

-949.8 

18.747 

-0.37017 


402 

-10.5 

0.209 

-0.00413 

55,914 

-1,161.7 

22.937 

-0.45290 


575 

-15.1 

0.299 

-0.00591 

58,951 

-1,269.9 

25.077 

-0.49516 


662 

-17.3 

0.345 

-0.00680 


( < P t R)‘ = 
(<5»r) 3 

IX - AqmjJ 


(66) 


Pitching-moment-distribution equations . - In the general case, the 
coefficients in the equation for pitching moment due to distributed wing 
lift are obtained by integrating the shear coefficients due to distributed 
lift from wing tip to center line in the x-direction. However, when the 
wing-quarter-chord locus is a straight line, this integration is equiva- 
lent to multiplying the corresponding bending-moment coefficients by 

-tan A /» . 
c/4 


The pitching-moment-equation coefficients due to external-store load 
are calculated by premultiplying the external-store load-coefficient mat- 
rix by a diagonal matrix whose elements are the chordwise distances from 
the store center of pressure to the vortex stations inboard of the store. 
For the sample problem, this moment-arm matrix is 


l/^esj “ 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

141.7 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

93.0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

44.3 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

4.4 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

- 28.81 


(67) 
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The resulting equation for pitching moment per unit of wing-root angle 
of attack is as follows: 


Load on wing: Load on external store: 


- 237 

5.1 

-0.102 

0.00201 


0 

0 

0 

0 



-1,810 

39.1 

-0.776 

0.01532 


0 

0 

0 

0 



-5,894 

126.0 

-2.497 

0.04928 


0 

0 

0 

0 



-12,994 

272.4 

-5.392 

0.10659 


0 

0 

0 

0 



qmp 

-23,447 

478.1 

-9.449 

0.18659 


0 

0 

0 

0 



O ^) 2 

-37,505 

738.6 

-14.570 

0.28754 

+ 

0 

0 

0 

0 


- 

(i^r) 5 

-55,395 

1,046.7 

-20.606 

0.40655 


2,024 

-25.6 

0.486 

-0.00958 




-77,302 

1,39**. 2 

-27.391 

0.54006 


■1,328 

-16.8 

0.319 

-0.00629 



1 - 

-103,384 

1,772.1 

-34.7^9 

0.68507 


633 

-8.0 

0.152 

-0.00299 



-133,769 

2,172.1 

-42.519 

0.85818 


63 

•0.8 

0.015 

-0.00030 



-150,051 

2,331-3 

-46.488 

O.91659 


-411 

5.2 

-0.099 

0.00195 




The pitching moment due to inertia twist per unit of load factor is 
given by the following equation: 



Load on wing: Load on external store: 


-50 

2.0 

-0.039 

O.OOO76 


0 

c 

0 

0 



-542 

14.9 

-0.294 

0.00579 


0 

c 

0 

0 



-1,991 

47.9 

-0.945 

0.01864 


0 

c 

0 

0 



-4,578 

103.4 

-2.038 

0.04024 


0 

c 

0 

0 



® n R 

-8,260 

181.2 

-3-570 

0.07050 


0 

c 

0 

0 



(< 3^)2 

-12,893 

279.2 

-5-504 

0.10868 

+ 

0 

c 

0 

0 



(qn ^)3 • 

-18,296 

394.7 

-7.783 

0.15369 


353 

-9.2 

0.184 

-0.00562 



(^r ) 4 

-24,319 

524.3 

-10.345 

0.20427 


252 

-6.1 

0.121 

-0.00258 



l - 

-30,842 

664 . 9 

-13-123 

0.25912 


no 

-2.9 

0.057 

-0.00115 



-37,740 

813.2 

-16.056 

0.31703 


11 

-0.3 

0.006 

-0.00011 



-41,266 

888.9 

-17.554 

0.34661 


-72 

1.9 

-0.037 

0.0007^ 

! 



Evaluation of - Equations have been written for shear, 

1 - Aqm R 


bending moment, and pitching moment in terms of qm R and . 

1 - Aqm^ 

In the evaluation of — 

1 


— given in equation ( 23 ), the quantity A 


is seen to be 
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c e,N+i 

c e,N 


( 70 ) 


As recommended in the "Method" section, values of C„ „ , and C a „ 

0 ,N +1 0 ,N 

corresponding to the 65 -percent-semispan station are used. From table IV 
which is used for the calculation of the C e j- coefficients, it can be 
seen that 


and 


C 6 ,N +1 



0.000000170540 


c e,N = { c 0,3f = -0.0000086353 

from which 


A = -0.0197306 


Additional Airloads 

The calculations just completed may now be used to determine the 
total additional airload shear, the bending moment and pitching moment, 
the horizontal-tail load per unit of normal load factor, the wing center 
of pressure, and the aerodynamic-center position. The numerical rela- 
tionship must be established between the wing-root angle of attack Arr 

and the airplane load factor An by use of equation (31) • The first 
s tep in the procedure is the calculation of wing-root loading components 
for selected values of the parameter qmp. 

Calculation of wing-root loading components .- The values of r , 

M b,a,r> M Y,a,r> L n,r> M b,n,r> and M Y,n,r (wing-root values) to be 
used in the succeeding calculations are evaluated from equations (62), 
(65)7 (68), (63), (66), and (69), respectively, by using only the bottom 
row of coefficients from each equation. These quantities, with the 
value of qmpj for which they are calculated, are presented in table VI. 

Additional wing and tail airloads and wing centers of pressure . - The 
quantities listed in table VI are used in the sample calculations pre- 
sented in table VII. The first seven columns of table VII represent the 

Acl 

steps involved in calculating — £ by use of equation (31) . Column 8 

An 
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lists the resulting values of Column 9 represents the evaluation 

of wing-root shear per unit of load factor from equation (33) • Column 10 
lists the horizontal-tail load per unit of load factor from equation ( 32 )* 

The wing-root pitching moment and bending moment per unit of load 
factor are listed in table VII in columns 13 and 14, respectively, as 
evaluated from equations of the same form as that of equation ( 33 )* The 
chordwise and spanwise centers of pressure x C p and y C p are evaluated 

in columns 15 and l 6 , respectively, by dividing the proper moment per 
unit of load factor by the shear per unit of load factor. Although the 
centers of pressure at qm^ = 0 appear indeterminate, they are equiva- 
lent to the centers of pressure for the theoretical rigid wing which may 
be calculated from the first coefficients in the rows of equations ( 62 ), 
( 65 ), and ( 68 ) corresponding to the wing-root position for both wing and 
external-store parts. For the sample problem. 


and 



- 1 ^ 1 ~ 411 = - 206.6 
714.0 + 14.3 



214359 + 3798 _ pqq = 
714.0 + 14.3 990 


Column 17 shows the incremental change in chordwise center of pressure 
from that at qmp = 0 and is expressed as a percentage of the wing mean 

aerodynamic chord. Thus, column 17 represents the change in wing- 
aerodynamic-center location due to flexibility. 

Although the actual location of the aerodynamic center of the air- 
plane may be different from that of the wing because of fuselage moments 
which were unaccounted for in the analysis, the change in aerodynamic- 
center location with change in qm^ should be practically unaffected 

by fuselage moments. If the aerodynamic -center location on a rigid wind- 
tunnel model is known, that location may he taken to be the actual value 
at qm^ = 0. The aerodynamic -center local ion at any other value of qm^ 

can be determined with good accuracy by acding the calculated change in 
aerodynamic-center location to the measured rigid-airplane position. 

Although only the shears and bending moments per unit of load factor 
at the center-line station were listed, those at other wing stations 
might have been calculated to obtain distributions along the semi span 
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and variations of the distribution with qm^. Figure 3(a) shows the 

shear distributions along the semispan per unit of load factor at values 
of qmR of 10 and 50 lb/ft 2 -deg. Each distribution was divided by the 

center-line shear. These ratios were plotted to show more clearly the 
fraction of the total aerodynamic lift carried outboard of any wing sta- 
tion as it varies with qrn^. The corresponding bending moments per unit 
of shear at the center line are plotted in figure 3(b). 

A plot such as that shown in figure 4 is useful in applying the 
calculations to design conditions of dynamic pressure and Mach number. 
The parameter qmp is plotted against q for several constant Mach 

numbers. The values of m^ used in this plot were obtained in wind- 
tunnel tests of the particular airplane used as an example. 


DISCUSSION 


The matrix equation expressing the change in angle -of -attack dis- 
tribution in terms of the applied angle-of -attack distribution is well 
suited to solution by matrix iteration. The sample problem included in 
this report illustrates that only three or four premultiplications of 
a column by a square matrix will give the elements of a simple expres- 
sion for final angle of attack resulting from any given initial angle 
of attack. This expression is readily evaluated for any value of the 
parameter qm^. The final spanwise distributions of lift, shear, 

bending moment, and pitching moment are also easily calculated. The 
solution can be carried to any desired degree of accuracy, and the 
degree of accuracy can be easily checked. 

Combining the iterative solution with the use of the downwash- 
f actor tables to obtain the matrix of aerodynamic -influence coefficients 
allowed the results of the sample problem to be calculated in about one- 
tenth of the estimated time required to obtain the solutions by calcu- 
lating the elements of pE>J and by using matrix inversion. 


The conditions under which the method of this report may be used 

,-l I 


are that the inverse of 




“r 


must exist and that the 


dominant latent root of this matrix must be real and distinct from its 
other latent roots. These conditions will be met in any practical prob- 
lem in aeroelasticity. The convergence of an infinite series of 


{ c e,Nf ( q - n ^) N 


terms is not a necessary condition, as should be evident 
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from the derivation of the equations pertinent to the iteration procedure 
and f ran the fact that solutions were obtained in the sample problem for 
values of qm^ above the value at which divergence of an infinite series 

would begin ^qmp = 50 lb/ft^-deg) . Actual wing divergence is indicated 

if the value of i - becomes infinite at some positive value of qm~ 

1 - Aqm^ 

This condition will occur only when the latent root A is positive, as 
it would be for a sweptforward wing. 

The methods of this report are designed to save time and effort in 
solving steady-state aeroelastic problems without degrading the accuracy 
of the results. There is, however, an added flexibility in that the 
various quantities are obtained as functions of the parameter qm^ and 

are readily applied to any flight condition of dynamic pressure and Mach 
number . This flexibility may make the iterative solution desirable even 
when automatic computing machines are available for matrix inversion. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., January 30, 1959- 
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APPENDIX A 

AERODYNAMIC -INFLUENCE COEFFICIENTS FROM 
TABLES OF DOWNWASH FACTOR 


The derivation of a system of aerodynamic -influence coefficients 
has been presented in reference 1. The computations involved in deter- 
mining any one element of the aerodynamic-influence-coefficient matrix 
[SiJ are rather lengthy, especially for manually operated desk calcu- 
lators. Tables and charts for obtaining downwash factors, based on 
lifting-line theory, have been in existence for some time in reference 2. 
An element of the [^l] matrix is obtainable from these tables by 

finding the downwash factor at the three-quarter-chord station in ques- 
tion and multiplying this factor by the reciprocal of the semispan of 
the horseshoe vortex which causes the downwash. The resulting element 
is identical, within the accuracy of the tables, with that which would 
be calculated by the method of reference 1. 


The tables of downwash factor in reference 2 have been included as 
tables I and II of this report with certain changes; the "corrector" 
vortex columns have been replaced by interpolating factors designated 
as K. Also, the values of downwash factor included in this report have 
been recalculated and are believed to be correct. The intervals used 
in the tables are such that linear interpolation will preserve the 

accuracy to four decimal places. A column for ^ = 38 has been 

h 

added in table I so that 10 equally spaced vortex stations on each 
semispan may be used. 

The arguments for tables I and II are ^ and ^ which are. 

h h * 

respectively, the streamwise and spanwise distances from any vortex 
station to any downwash station, nondimens ionalized with respect to the 

Ax 

vortex semispan. The argument — is negative when the downwash sta- 
tion is ahead of the vortex station. 


It will be a help in using the tables if the values of 


^ and ^ 
h h 

are arranged in matrix form. A systematic method for determining these 
and |^~j matrices for 10 vortex semispan stations is given as 

follows; With the use of a coordinate system having the origin at the 
point of intersection of the wing quarter-chord line with the airplane 
center line, the nondimens ional lateral coordinate of any vortex center 
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is given by 



(Al) 


where N y is the number of vortices on both wings and Ti y is the vortex 

spanwise station expressed as a fraction of the semispan. For example, 
at station 1 in the sample problem, T) y = 0.95> Ky = 20 , and ^Q y = ^9- 


The coordinates of all the vortex stations expressed as a row matrix 
fZ\ I are, for 10 equally spaced vortex stations on each semispan, 

AWvJ 

= [19 17 15 15 11 9 7 5 5 lj 

I©! (ft (ft ••• (ftoj (A2> 


The nondimens ional lateral coordinates of the downwash control points 
are the same as those of the vortex centers. 


Two 


¥ 


matrices will be needed. The first 


i 

y 


Ajrl 


hJS 


, will con- 


tain the nondimens ional spanwise distances from each vortex center to 
every downwash control point on the same ving. The second, 1 will 

contain the nondimensional spanwise distances from each vortex station to 
every downwash control point on the opposite wing. The first row of the 


matrix 


Ay 


may be obtained by subtracting ^ of downwash control 

s 


For 10 equally 

Ay 
h 


point 1 from each element of the row matrix j^jQ v J 
spaced vortex stations on each semispan, *;he first row of 
be 

[19 17 15 13 ll 9 7 5 3 lj - [_19 19 19 19 19 19 19 19 19 19 J 

= [o 2 4 6 8 10 12 14 16 i8j 


will 


(A3) 


Absolute values are used since the downwash factor is independent of 
the sign of 
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The first row of 
obtain 


Ay 

_ h Jo 


is obtained by adding the two rows to 


[58 36 34 32 30 28 26 24 22 20j 


It will be seen that the matrices 


^ and ^1 can be obtained by 

LhJs LkJo 


finding the difference and the sum, respectively, of the two matrices 


U 


and 


fca 


19 

17 

15 

13 

ll 

9 

7 

5 

3 


1 

19 

17 

15 

13 

11 

9 

7 

5 

3 


1 

19 

17 

15 

13 

ll 

9 

7 

5 

3 


1 

19 

17 

15 

13 

ll 

9 

7 

5 

3 


1 

19 

17 

15 

13 

11 

9 

7 

5 

3 


1 

19 

17 

15 

13 

ll 

9 

7 

5 

3 


1 

19 

17 

15 

13 

11 

9 

7 

5 

3 


1 

19 

17 

15 

13 

11 

9 

7 

5 

3 


1 

19 

17 

15 

13 

11 

9 

7 

5 

5 


1 

19 

17 

15 

13 

ll 

9 

7 

5 

3 
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"19 

19 

19 

19 

19 

19 

19 

19 

19 

19” 


17 

17 

17 

17 

17 

17 

17 

17 

17 

17 


15 

15 

15 

15 

15 

15 

15 

15 

15 

15 


13 

13 

13 

13 

13 

13 

13 

13 

13 

15 


11 

ll 

11 

ll 

ll 

ll 

ll 

ll 

ll 

ll 


9 

9 

9 

9 

9 

9 

9 

9 

9 

9 


7 

7 

7 

7 

7 

7 

7 

7 

7 

7 


5 

5 

5 

5 

5 

5 

5 

5 

5 

5 


3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


l 

l 

l 

l 

1 

l 

1 

l 

1 

l 



(A4) 


(A5) 
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where 


[( 5 ) 


consists of 10 rows identical to the first row and 


[V- / V J 

consists of 10 columns identical to the first column. Thus, 



= 

f?r 



_ h _ 

S 

\ h /V 


\ E /d 


(a6) 


which can be expressed as 


and 


Ay 

hJs 


0 

2 

4 

6 

8 

10 

12 

14 

16 

18 


2 

0 

2 

4 

6 

8 

10 

12 

14 

16 


4 

2 

0 

2 

4 

6 

8 

10 

12 

14 


6 

4 

2 

0 

2 

4 

6 

8 

10 

12 


8 

6 

4 

2 

0 

2 

4 

6 

8 

10 


10 

8 

6 

4 

2 

0 

2 

4 

6 

8 


12 

10 

8 

6 

4 

2 

0 

2 

4 

6 


14 

12 

10 

8 

6 

4 

2 

0 

2 

4 


16 

14 

12 

10 

8 

6 

4 

2 

0 

2 


18 

l6 

14 

12 

10 

8 

6 

4 

2 

0 
















0 


1 ji 

& 




(A7) 


which can be expressed as 
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Ay 

h 


0 


38 36 

36 34 


34 32 30 28 26 24 22 


20 


32 30 28 26 24 22 20 18 


34 32 30 28 26 24 22 


20 


18 16 


32 

30 


30 28 26 24 22 


20 


18 16 14 


28 26 24 22 


20 


18 16 14 12 


28 26 24 22 


20 


18 16 14 12 


26 24 

24 22 


22 20 18 16 14 12 


20 


18 16 14 12 


22 20 18 16 14 12 


20 


18 16 14 12 


10 


10 

8 


10 

8 

6 


10 

8 

6 

4 


10 

8 

6 

4 

2 


If using fewer than 10 equally spaced vortex stations on each semi- 
span is desired, the matrices j and (fijd can mo< ^^* e( ^ ^o 

suit by omitting the required number of rows and columns from the top 
and left. 


By using the same coordinate origin as before, a matrix of the non- 

(h) v ] uhen 


dimensional longitudinal coordinates of the vortex centers 
the quarter-chord locus is a straight line is 




tan A 


c/4 


(A8) 


For the downwash control points on the three -quarter -chord line, 


l)c 


'flY 

\h) d 


where each element in the matrix 
divided by the vortex semispan h. 


tan 


c/4 + 


2h 


(A9) 


2h 


is half the wing-root chord 


The correct matrix of nondimens ional longitudinal distances between 


vortex and control points 


Ax 


is given by the equation 
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(A10) 


Ax _ 
h 


= (A 

w< 


tan a 5 c A + i - [(h) J tan A <=A 


(All) 


Only one y — J matrix is required since chordwise distances from 

a vortex to any control point are the same whether the vortex and control 
point are on the same wing or on opposite wings. 

Two complete matrices of downwash factors will he obtained: the 

first by using ^ and j^j and the second by using ^ ^ and 

as arguments in tables I and II. The symmetrical-load matrix 

is calculated by adding these two matrices of downwash factors and multi- 
plying each element by the reciprocal of the vortex semispan i. For 

antisymmetric loadings, the value of [sj is the difference between 
the two downwash-f actor matrices multiplied by 


Tables I and II are useful only for equally spaced horseshoe vor- 

Ay Ax 

tices. Unequal spacing would require interpolation in both and — , 

which is not practical. 

The columns in tables l(a) and II entitled "K" contain interpolating 
factors which were calculated by dividing each increment in downwash 

factor by the corresponding increment in Each factor is listed in 

the same row with the smaller of the two values used to calculate 

h 

it. In order to obtain the downwash facte r F corresponding to any 

value of — , find the value of downwash factor F' corresponding to the 
h 

next smaller value of in the table (in absolute value). Add to 

F ' the product of ^ and the interpolai ing factor and subtract the 

A Y * 

product of - — and the interpolating factor. Thus, 
h 


F = F' +^K -^1 
h n 


(A12) 
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As an example, suppose it is required that a value of F be 
obtained that corresponds to ^ = 1.3421 and ^ = 0. By using 

Ay Ay I 

table II, the next smaller value of =- is I .34 which equals =— for 

h h 

which the value of F ! is 4.4956 and the interpolating factor is -0.660. 
The required value of F is 

F = 4.4956 + 1.342l(-0.660) - 1,3400( -0.660) 

F = 4.4942 

With the inclusion of the interpolating factors, all information 
necessary to obtain the downwash factors corresponding to a given value 

/\y 

of — lies in one row of the tables, the row corresponding to the 

A y" 

value of that is next smaller than the given value. Interpolation 

h 

can be done without writing or remembering any intermediate steps with 
any calculating machine capable of accumulative multiplication. Where 
the interpolating factors have been omitted > 10 j , interpolation by 

inspection will be sufficiently accurate. 
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TABLE I.- DOWNWASH FACTOR F‘ 


(a) Values of F far £%/ h = 0 to ±10 


Ax/h 


= 0 


= +2 

ty/h 

= ±4 

/ty/h 

* ±6 

fy/h = ±8 

Zty/h 

= +10 

F 

1 K 

F 

K 

F 

K 

F 

K 

F 

K 

F 

K 

-30 



-0.0011 


-0.0011 


-o.oon 


-0.0011 


-0.0010 









-25 

-.0016 

-.0016 

-.0015 

-.0015 

-.0014 








-20 



-0.0025 

-0.00018 

-0.0024 

-0.00018 

-0.0023 

-0.00018 

-0.0022 


-0.0021 


-18 



-.0051 

-.00030 

-.0050 

-.00030 

-.0029 

-.00030 

-.0027 

-0 .00025 

-.0025 

-0.00020 

-16 



-.0059 

-.00040 

-.0037 

-.00040 

-.0035 

-.00030 

-.0053 

-.00030 

-.0030 

-.00025 

-14 



-.0050 

-.00055 

-.0048 

-.00055 

-.0045 

-.00050 

- .0041 

-.00040 

-.0037 

-.00035 

-12 



- .0068 

- .00090 

-.0064 

-.00080 

-.0059 

-.00070 

-.0053 

-.00060 

-.0046 

-.00040 

-10 



-0.0097 

-0.00145 

-0.0089 

-0.00125 

-0.0079 

-0 .00100 

-0 .0069 

-0.00080 

-0.0059 

-0 .00060 

-9 



-.0119 

-.00220 

-.0108 

-.00190 

-.0093 

-.00140 

-.0079 

-.00100 

-.0066 

-.00070 

-8 



- .0149 

- .00500 

-.0132 

-.00240 

-.0111 

-.00180 

- .0092 

-.00130 

-.0075 

-.00090 

-7 



- .0192 

-.00450 

- .0165 

-.00530 

- .0134 

-.00230 

- .0107 

-.001 l ;Q 

-.0086 

- .00110 

-6 



-.0256 

- .00640 

- .0210 

-.00450 

- .0163 

-.00290 

- .0125 

-.00180 

-.0097 

-.00110 

-5 



-0.0356 

-0.01000 

-0.0274 

-0.00640 

-0.0201 

-0.00380 

-0.0148 

-0.00230 

-0.0111 

-0.00140 

-4.5 



-.0429 

- .01460 

-.0317 

-.00860 

- .0224 

-.00460 

-.0160 

- .00240 

-.0119 

-.00160 

-4.0 





-.0368 

-.01020 

- .0249 

-.00500 

- .0174 

- .00280 

-.0126 

- .00140 

-3.5 





-.0431 

-.01260 

-.0278 

-.00580 

-.0189 

- .00300 

-.0135 

-.00180 

-5-0 





-.0507 

-.01520 

-.0311 

-.00660 

-.0205 

- .00320 

-.0143 

-.00160 

-2.5 





- .0599 

-.01841 

- .0347 

-.00720 

- .0222 

- .00340 

--0153 

- .00200 

-2.0 

(2) 



-0.0709 

-0.02200 

-0.0386 

-0.00780 

-0.0239 

-0.00340 

-0 .0162 

-0.00180 

-1.8 





-.0759 

-.02500 

- .0403 

-.00850 

- .0247 

-.00400 

-.01 66 

-.00200 

-1.6 





-.0812 

-.02650 

- .0420 

-.00850 

-.0254 

-.00350 

-.0170 

-.00200 

-1.4 





-.0869 

-.02850 

- .0438 

-.00900 

- .0262 

-.00400 

- .0174 

-.00200 

-1.2 





-.0928 

-.02950 

-.0456 

- .00900 

-.0270 

- .00400 

-.0178 

-.00200 

-1.0 





- .0990 

-.03100 

-.0475 

-.00950 

-.0278 

-.00400 

-.0102 

-.00200 

-.8 





- . 1056 

- .03300 

- .0494 

-.00950 

-.0285 

-.00350 

-.0186 

-.00200 

-0.6 





-0.1123 

-0.03350 

-0.0513 

-0.00950 

-0.0293 

-0.00400 

-0.0190 

-0.00200 

-.5 





- . 1156 

-.03300 

-.0523 

- .01000 

-.0298 

-.00500 

-.0192 

-.00200 

-.4 





- . 1192 

-.03600 

-.0532 

- .00900 

-.0502 

- .00400 

-.0194 

-.00200 

-0.5 





-0 . 1227 

-0.03500 

-0.0542 

-0.01000 

-0.0305 

-0.00300 

-0.0196 

-0.00200 

-.25 





- .1245 

-.03600 

-.0547 

-.01000 

-.0308 

-.00600 

-.0197 

-.00200 

- .20 





-.1262 

- .03400 

- .0552 

- .01000 

-.0310 

-.00400 

-.0198 

-.00200 

-.15 





-.1280 

-.03600 

- .0557 

-.01000 

-.0312 

-.00400 

-.0199 

-.00200 

-.10 





-.1298 

-.03600 

- .0561 

-.00800 

-.0314 

-.00400 

-.0200 

-.00200 

-.05 





-.1315 

-.03400 

- .0565 

- .00800 

-.0316 

-.00400 

-.0201 

-.00200 

0 




-0.1333 

-0 .03600 

-0.0571 

-0.01200 

-0.0318 

-0.00400 

-0.0202 

-0.00200 

0.05 





-0.1351 

-0.05400 

-0.0577 

-0.00800 

-0.0320 

-0.00400 

-0.0203 

-0.00200 

.10 





-.1368 

-.03800 

- .0581 

- .01000 

-.0322 

-.00400 

-.0204 

-.00200 

.15 





-.1387 

-.03400 

-.0586 

-.01000 

-.0324 

-.00400 

-.0205 

- .00200 

.20 





-.1404 

-.03600 

-.0591 

-.01000 

-.0326 

-.00400 

-.0206 

-.00200 

.25 





-.1422 

- .03600 

-.0596 

- .01000 

-.0328 

-.00400 

-.0207 

-.00200 

.5 





-.1440 

-.03500 

-.0601 

- .00900 

-.0530 

-.00400 

-.0208 

-.00200 

0.4 





-0.1475 

1 -0.03500 

-0 .0610 

-0.01000 

-0.0334 

-0.00400 

-0.0210 

1 

-0.00200 

-5 





-.1510 

-.03400 

-.0620 

- .01000 

-.0336 

- .00400 

-.0212 

-.00200 

.6 





-.1544 

-.03350 

-.0630 

-.00950 

-.0342 

- .00400 

-.0214 

-.00200 

0.8 





-0.l6ll 

-0 .03250 

-0.0649 

-0.00950 

-0.0350 

-0.00350 

-0.0218 

-0.00200 

1.0 




- . 1676 

-.03150 

-.0668 

-.00950 

-.0357 

- .00400 

-.0222 

-.00200 

1.2 





-.1739 

-.02950 

-.0687 

-.00900 

-.0365 

-.00400 

- .0226 

-.00200 

1.4 





-.1798 

-.02800 

-.0705 

- .00850 

-.0373 

-.00400 

-.0230 

-.00200 

1.6 





-,1854 

- .02700 

- .0722 

-.00900 

-.0381 

-.00350 

- .0234 

-.00200 

1.8 





-.1908 

- .02450 

-.0740 

-.00850 

-.0388 

-.00400 

-.0238 

-.00200 

2.0 





-.1957 

- .02220 

-.0757 

- .00780 

-.0396 

-.00360 

-.0242 

-.00180 

2.5 





-0.2068 

-0.01840 

-0.0796 

-0.00720 

-0.0414 

-0.00320 

-0.0251 

-0.00200 

3-0 





-.2160 

-.01520 

- .0832 

-.00660 

-.0430 

- .00320 

- .0261 

-.00160 

3-5 





-.2236 

- .01240 

-.0865 

-.00580 

-.0446 

-.00300 

.-.0269 

-.00180 

4.0 





- .2298 | 

- .01040 

-.0894 

-.00500 

-.0461 

- .00280 

-.0278 

-.00160 

4.5 

4.0488 

0.01840 

-1.2905 

-0.01440 

- .2350 j 

-.00840 

-.0919 

-.00460 

-■0475 

-.00240 

-.0286 

-.00140 

5 

4.0596 

.01200 

-1.2977 

- .01010 

-.2392 

-.00650 

-.0942 

-.00380 

- .0487 

-.00230 

-.0293 

-.00140 

6 

4.0276 

0.00730 

-1.5078 

-0 .00640 

-0.2457 

-0.00450 

-0.0980 

-0.00290 

-0.0510 

-0.00180 

-0 .0307 

-0.00120 

7 

4.0205 

.00470 

-1.5142 

-.00450 

-.2502 

-.00350 

-.1009 

-.00230 

- .0528 

-.00150 

-.0319 

-.00100 

8 

4.0156 

.00550 

-1.5185 

-.00500 

- .2535 

-.00240 

-.1032 

-.00180 

-.0543 

-.00130 

-.0329 

-.00090 

9 

4.0125 

.00250 

-1.5215 

-.00220 

- .2559 

-.00180 

-.1050 

-.00140 

-.0556 

-.00100 

- .0338 

-.00070 

10 

4.0100 

.00155 

-1.3237 

-.00140 

- .2577 

- .00130 

-.1064 

-.00100 

-.0566 

-.00085 

-.0345 

-.00065 

12 

4.0069 

0.00095 

-1.5265 

-0.00090 

-0.2603 

-0.00080 

-0.1084 

-0.00070 

-0.0583 

-0.00050 

-0.0358 

-0.00045 

14 

4.0050 

.00055 

-1.5285 

-.00060 

-.2619 

-.00050 

-.1098 

-.00050 

-.0594 

-.00040 

- .0367 

-.00035 

16 

4.0039 

.00040 

-1.3295 

-.00040 

-.2629 

-.00040 

-.1108 

-.00030 

-.0602 

-.00030 

- .0374 

-.00025 

18 

4.0051 

.00050 

-1.3305 

-.00050 

-.2637 

-.00025 

-.1114 

-.00025 

-.0608 

-.00025 

-.0379 

-.00020 

20 

4.0025 

.00018 

-1.3309 

-.00018 

-.2642 

-.00018 

-.1119 

-.00010 

- .0613 


- .0383 


25 

4.0016 


-1.5318 


-0.2651 


-0.1128 


-0.0620 


-0.0390 


30 







4.0011 

-1.5322 

-.2656 

-.1132 

-.0624 

-.0394 








■’■Values of F are obtained from reference 2 with slight modification. 
"Values in this range are not used. 

^Values in this range are given in table II. 




TABLE I.- DOWNWASH FACTOH F - Coi eluded 


(b) Values of F for ty/h = +12 to +38 








F for values of <ty/h 

of - 





! 


+ 12 

+ 14 

+16 

+18 

+ 20 

±22 

±24 

+26 

±28 

±30 

±32 

±34 

136 

±58 

-30 

-0.0010 

-0.0010 

-0.0009 

-0.0009 

-0.0008 

-0.0008 

-0.0008 

-0.0007 

-0.0007 

-0.0007 

-0.0006 

-0.0006 

-0.0006 

-0.0006 

-^5 

-.0014 

-.0013 

-.0012 

-.0012 

-.0011 

-.0010 

-.0010 

-.0009 

-.0009 

-.0008 

-.0007 

-.0007 

-.0007 

-.0007 

-i^O 

-0.0020 

-0.0019 

-0.0017 

-0.0016 

-0.0013 

-0.0013 

-0.0012 

-0.0012 

-0.0011 

-0.0010 

-0.0009 

-0.0009 

-0.0008 

-0.0008 

-10 

-.0023 

-.0022 

-.0020 

-.0018 

-.0017 

-.0015 

-.0014 

-.0013 

-.0012 

-.0011 

-.0010 

-.0009 

-.0009 

-.0008 

-lb 

-.0028 

-.0029 

-.0023 

-.0021 

-.0019 

-.0017 

-.0015 

-.0014 

-.0013 

-.0012 

- .0011 

-.0010 

-.0009 

-.0008 

-14 

-.0034 

-.0030 

-.0027 

-.0024 

-.0021 

-.0019 

-.0017 

-.0016 

-.0014 

- .0013 

-.0012 

-.0011 

-.0010 

-.0009 

-12 

-.0041 

-.0036 

- .0031 

-.0028 

-.0024 

-.0022 

-.0019 

-.0017 

- .0015 

- .0014 

- .0013 

-.0012 

-.0011 

-.0010 

-10 

-0.0090 

-0.0043 

-0.0037 

-0.0032 

-0.0028 

-0.0024 

-0.0021 

-0.0019 

-0.0017 

-0.0015 

-0.0014 

-0.0012 

-0.0011 

-0.0010 

-9 

-.0096 

-.0047 

-.0040 

-.0034 

-.0030 

-.0026 

-.0023 

-.0020 

-.0018 

-.0016 

-.0014 

-.0013 

-.0012 

- .0010 

-8 

-.0062 

- .0092 

- .0043 

-.0037 

-.0051 

-.0027 

-.0024 

-.0021 

-.0019 

- .0017 

-.0015 

-.0013 

-.0012 

- .0012 

-7 

-.0069 

- .0097 

-.0047 

- .0039 

-.0034 

-.0029 

-.0025 

-.0022 

- .0019 

-.0017 

-.0015 

- .0014 

- .0012 

- .0012 

-- 

-.0077 

-.0062 

-.0051 

-.0042 

-.0036 

-.0051 

-.0026 

-.0023 

-.0020 

-.0018 

-.0016 

-.0014 

-.0013 

- .0012 


-0.0086 

-0.0068 

-0.0053 

-0.0045 

-0.0038 

-0.0032 

-0.0028 

-0.0024 

-0 .0021 

-0.0019 

-0.0017 

-0.0015 

-0.0013 

-0.0012 

-4.5 

-.0091 

-.0071 

-.0057 

- .0047 

-.0039 

-.0033 

-.0028 

- .0025 

-.0021 

-.0019 

-.0017 

-.0015 

-.0014 

-.0012 

-4.0 

-.0099 

-.0074 

-.0039 

-.0048 

-.0040 

-.0034 

-.0029 

-.0025 

- .0022 

-.0019 

-.0017 

-.0015 

-.0014 

-.0013 

-.5.9 

-.0101 

- .0078 

-.0062 

-.0030 

-.0041 

-.0053 

- .0030 

-.0026 

- .0022 

-.0020 

-.0017 

-.0016 

- .0014 

-.0013 

-5.0 

-.0106 

-.0081 

-.00? .4 

-.0052 

-.0043 

-.0036 

-.0030 

- .0026 

- .0023 

- .0020 

-.0018 

-.0016 

-.0014 

-.0014 

-2.9 

- .0111 

- .0063 

-.00f*.. 

-.0053 

-.0044 

-.0037 

-.0031 

- .0027 

-.0023 

-.0020 

-.0018 

-.0016 

-.0014 

-.0014 

-2.0 

-0.0117 

-0.0088 

-0.0069 

-0.0055 

-0.0045 

-0.0038 

-0.0032 

-0.0027 

-0.0024 

-0 .0021 

-0.0010 

-0.0016 

-0.0015 

-0.0014 

-1.8 

- .0119 

-.0089 

- .0070 

-.0056 

- .0046 

-.0038 

-.0032 

- .0028 

-.0024 

-.0021 

-.0018 

-.0016 

-.0015 

-.0014 

-1.6 

- .0121 

- .0091 

- .0071 

- .00% 

- .0046 

-.0038 

- .0032 

- .0028 

-.0024 

-.0021 

-.0019 

-.0017 

-.0015 

-.0014 

- ! .4 

-.0124 

-.0092 

- .0072 

-.0057 

- .0047 

-.0039 

-.0033 

- .0028 

-.0024 

-.0021 

-.0019 

-.0017 

-.0015 

-.0014 

-1 .2 

- .0126 

-.0094 

-.0073 

-.0090 

- .0047 

-.0039 

-.0033 

-.0028 

- .0024 

-.0021 

-.0019 

-.0017 

-.0015 

-.0014 

-1 .0 

-.0128 

- .0099 

-.0074 

-.0058 

- .0048 

-.0040 

-.0033 

- .0028 

-.0025 

- .0022 

-.0019 

-.0017 

-.0015 

-.0014 

-.8 

- .0131 

- .0097 

- .0075 

-.0059 

-.0048 

-.0040 

-.0034 

- .0029 

-.0025 

-.0022 

-.0019 

-.0017 

- .0015 

-.0014 

-O.t, 

-0.0135 

-0.0098 

-0.0075 

-0.0060 

-0.0049 

-0.0040 

-0 .0034 

-0.0029 

-0.0025 

-0.0022 

-0.0019 

-0.0017 

-0.0015 

-0.0014 

- .5 

- .01 34 

-.0099 

-.0076 

-.0060 

-.0049 

- .0040 

-.0034 

-.0029 

-.0025 

-.0022 

-.0019 

-.0017 

-.0015 

-.0014 

- .1* 

-.0139 

-.0100 

-.0076 

-.0061 

-.0049 

- .0041 

-.0034 

-.0029 

-.0025 

- .0022 

-.0019 

-.0017 

-.0015 

-.0014 

-0.3 

-0.0130 

-0.0100 

-0.0077 

-0.0061 

-0.0049 

-0.004 j. 

-0.0034 

-O.O029 

-0.0025 

-0.0022 

-0.0019 

-0.0017 

-0.0015 

-0 .0014 

- .29 

-.0157 

-.0101 

-.0077 

-.0061 

-.0050 

-.0041 

-.0034 

-.0029 

-.0025 

- .0022 

-.0019 

- .0017 

-.0015 

-.0014 

- ...O 

-.0157 

-.0101 

-.0077 

-.0061 

-.0050 

-.0041 

-.0054 

-.0029 

-.0025 

- .0022 

-.0019 

-.0017 

- .0015 

- .0014 

- .19 

- .0158 

-.0102 

-.0076 

-.0061 

- .0050 

-.0041 

-.0035 

-.0029 

-.0025 

- .0022 

-.0019 

-.0017 

- .0015 

- .0014 

- .10 

- .0139 

-.0102 

-.0078 

- .0062 

- .0050 

-.0041 

-.0035 

-.0030 

-.0025 

-.0022 

-.0019 

-.0017 

-.0015 

- .0014 

- .09 

-.0140 

-.0102 

- .0078 

-.0062 

-.0050 

-.0041 

-.0035 

-.0030 

-.0025 

- .0022 

- .0020 

-.0017 

-.0015 

_| 

-.0014 

0 

-0.0140 

-0.0103 

-0.0078 

-0.0062 

-0.0050 

-0.0041 

-0.0035 

-0.0030 

-0.0026 

-0 .0022 

-0.0020 

-0 .0017 

1 

-0.0015 

-0 .0014 

0.09 

-0.0140 

-0.0103 

-0 .0079 

-0 .0062 

-0.0050 

-0 .0042 

-0.0035 

-0.0030 

-0.0026 

-0.0022 

-0.0020 

-0.0017 

-0.0015 

-0-0014 

.10 

-.0141 

- .0103 1 

- .0079 

-.0062 

-.0050 

- .0042 

-.0035 

-.0030 

- .0026 

- .0022 

-.0020 

-.0017 

I -.0015 

-.0014 

.19 

-.0142 

- .0104 

-.0079 

-.0062 

-.0050 

-.0042 

-.0035 

-.0030 

-.0026 

- .0022 

- .0020 

-.0017 

- .0016 

-.0014 

.20 

-.0142 

-.0104 

-.0079 

-.0063 

- .0031 

-.0042 

-.0035 

-.0030 

-.0026 

- .0022 

-.0020 

-.0017 

-.0016 

-.0014 

.29 

-.0143 

- .0104 

-.0080 

-.0063 

-.0051 

-.0042 

-.0035 

-.0050 

- .0026 

- .0022 

-.0020 

-.0017 

- .0016 

-.0014 

■ 5 1 

- .0144 

- .0109 

- .0080 

-.0065 

-.0051 

-.0042 

-.0035 

- .0030 

-.0026 

- .0022 

- .00 20 

-.0017 

- .0016 | 

-.0014 

0.4 

-0.0149 

-0 .0109 

-0.0000 

-0.0063 

-0.0051 

-0.0042 

-0.0035 

-0.0030 

-0.0026 

-0.0023 

-0.0020 

-0.0018 

-0.0016 

-0.0014 

.9 

- .0146 

- .0106 

-.0081 

-.0064 

-.0031 

-.0042 

- .0036 

- .0030 

- .0026 

-.0023 

-.0020 

-.0018 

- .0016 

- .0014 

.6 

- .0147 

- .0107 

- .0081 i 

-.0064 

-.0052 

-.0043 

- .0036 

-.0030 

-.0026 

-.0023 

-.0020 

-.0018 

- .0016 

-.0014 

0.8 

-0 .0149 

-0.0109 

-0.0082 

-O.OO65 

-0 .0052 

-0.0043 

-0.0036 

-0.0031 

-0.0026 

-0.0023 

-0.0020 

-0 .0018 

-0-0016 

-0.0014 

1.0 

- .0192 

- .0110 

-.0003 

-.OO65 

-.0055 

- .0043 

- .0036 

- .0031 

1 - .0026 

-.0023 

-.0020 

-.0018 

-.0016 

-.0014 

1.2 

-.0154 

-.0111 

-.0084 

-.0066 

-.0053 

-.0044 

- .0037 

- .0031 

[ -.0027 

-.0023 

-.0020 

- .0018 

-.0016 

-.0014 

1.4 

- .0196 

-.0113 

-.0085 

-.0067 

-.0054 

-.0044 

-.0037 

-.0031 

1 - .0027 

-.0023 

-.0020 

- .0018 

-.0016 

-.0014 

1 .6 

-.0199 

-.0114 

-.0086 

-.0067 

-.0054 

- .0044 

-.0037 

-.0031 

- .0027 

- .0023 

-.0021 

-.0018 j 

-.0016 

-.0014 

1.8 

- .0161 

-.0116 

-.0087 

-.0068 

-.0053 

-.0045 

-.0037 

- .0032 

- .0027 

-.0024 

-.0021 

- .0018 

-.0016 

-.0014 

2.0 

-.0163 

-.0117 

-.0088 

-.0069 

-.0039 

-.0045 

- .0038 

-.0032 

- .0027 

-.0024 

-.0021 

- .0018 

-.0016 

-.0014 

2.9 

-0.0169 

-0.0121 

-0.0091 

-0.0070 

-0.0056 

-0.0046 

-0.0038 

-0 .0032 

-0.0028 

-0 .0024 

-0 .0021 

-0.0019 

-0-0017 

-0 .0014 

5.0 

- .0174 

- .0124 

-.0093 

-.0072 

-.0058 

-.0047 

-.0059 

-.0033 

- .0028 

-.0024 

-.0021 

-.0019 ! 

-.0017 

-.0014 

5.9 

- -0179 

-.0120 

-.0095 

- .0074 

-.0059 

- .0048 

-.0040 

-.0034 

-.0029 

-.0025 

-.0022 

-.0019 

-.0017 

-.0015 

4.0 

-.0184 

-.0131 

- .0098 

-.0075 

-.0060 

-.0049 

- .0041 

-.0034 

-.0029 

-.0025 

-.0022 

-.0019 

-.0017 

-.0015 

4 .9 

-.0189 

-.0134 

- .0100 

-.0077 

-.0061 

-.0050 

-.0041 

-.0035 

-.0030 

-.0026 

-.0022 

-.0020 

-.0017 

-.0016 | 

5 

-.0194 

- .0137 

- .0102 

-.0079 

-.0062 

-.0051 

-.0042 

-.0035 

-.0030 

-.0026 

-.0023 

-.0020 ' 

-.0018 

-.001.6 

o 

-0.0203 

-0 .014 5 

-0 .0106 

-0.0082 

-O.OO69 

-0.0052 

-0.0043 

-0.0036 

-0.0031 

-0.0027 

-0.0023 

-0.0020 

-0.0018 

-0.0016 

7 

- .0211 

-.0149 

- .0110 

-.0084 

-.0067 

-.0054 

-.0045 

-.0037 

-.0032 

-.0027 

-.0024 

-.0021 

-.0018 

-.0016 

6 

-.0210 

-.0194 

- .0114 

-.OO87 

-.0069 

-.0056 

-.0046 

-.0038 

-.0033 

-.0023 

-.0024 

-.0021 

-.0019 

-.0016 

9 

- .0224 

- .0138 

- .0117 

-.0090 

-.0071 

-.0057 

-.0047 

-.0039 

-.0033 

-.0029 

-.0025 

-.0022 

-.0019 

-.0017 

10 

-.0230 

- .0162 

-.0120 

-.0092 

-.0073 

-.0059 

-.0048 

-.0040 

- .0034 

-.0029 

-.0025 

-.0022 

-.0020 

-.0017 

12 

-0.0239 

-0 .0170 

-0.0126 

-0.0096 

-0.0076 

-0.006l 

-0.0050 

-0.0042 

-0.0036 

-0.0031 

-0.0026 

-0.0023 

-0.0020 

-0.0018 

14 

-.0246 

- -0179 

-.0130 

-.0100 

-.0079 

-.0064 

-.0052 

-.0044 

--0037 

-.0032 

-.0027 

-.0024 

-.0021 

-.0019 

16 

-.0292 

-.0180 

-.0134 

-.0103 

-.0081 

-.0066 

-.0054 

-.0045 

-.0030 

-.0033 

-.0028 

-.0025 

-.0022 

-.0020 

18 

-.0296 

-.0184 

-.0137 

-.0106 

-.0084 

-.0068 

-.0056 

-.0047 

--0039 

-.0034 

-.0029 

-.0025 

-.0022 

-.0020 

20 

-.0260 

- .0187 

- .0140 

- .0108 

-.0086 

-.0069 

-.0057 

-.0048 

- .0040 

-.0035 

-.0030 

-.0026 

-.0023 

-.0020 

29 

-0.0266 

-0.0192 

-0 .0145 

-0.0112 

-0.0009 

-0 .0073 

-O.OObO 

-0.0050 

-0.0043 

-0 .0037 

-0.0032 

-0.0028 

-0.0024 

-0.0021 

30 

-.0270 

-.0196 

-.0148 

- .0115 

-.0092 

-.0075 

-.0062 

-.0052 

: -.0044 

-.0038 

-.0033 

-.0029 

-.0025 

-.0022 



TABU; II.- DOWNWASH FACTOR F FOR / h = 0 and ±2 


j^KxpntiLi t on of i,«*i’tion:i (5) of tnhlt- l(ri)j 



Ay./*! = 

0 




= t? 

F 

K 

Ax/h 

F 

K 


Ax/h 

F 

K 

I2.L976 

11.7321 

-46.980 

-42.420 

0.85 

.86 

9 .0883 

-2.060 


-4.0 

-0 .0526 
- .0571 


9.0677 

-2.050 


-3-8 

-0 .0230 

11.3079 

-30.970 

.87 

9-0472 

-1.960 


-3.6 

-.0629 

- .0290 

10.9222 

-39.230 

.88 

9.0276 

-1.930 


-3-4 

-.0690 

-.0305 

10.^699 

-32.420 

.89 

9.0083 

-1.870 


-3.2 

- .0759 

-.0349 

10.2497 

-29.820 

■ 90 

4 .9896 

-1.800 


-3-0 

-.0839 

-.0400 

9-9479 

-27.460 

•91 

4-9716 

-1 .740 


-2.6 

-.0933 

-.0470 

9.6729 

-29.490 

-92 

4 .9942 

-1.740 


-2.6 

-.1042 

-.0945 

9-4180 

-23.780 

• 93 

4 .9368 

-1 .670 


-2.4 

-.1167 

-.0625 

9 . 1002 

-21.990 

• 94 

4 .92OI 

-1 .640 


-2.2 

-.1316 

-.0745 

8.9605 

-20.610 

.95 

4-9037 

-1.570 


-2.0 

-.1494 

-.0890 

8.7942 

-19-180 

■96 

4.0880 

-1.560 


-1.9 

-.1595 

-.1010 

8.9624 

-18.010 

•97 

4.8724 

-1.490 


-1.8 

- . 1706 

-.1110 

8.3823 

-16.870 

■98 

4.8975 

-1.500 


-1.7 

-.1825 

-.1190 

8.2136 

-19-980 

•99 

4.8429 

-1 .400 


-1.6 

-.1956 

-.1310 

8.0938 

-14 .920 

1.00 

4 .8205 

-1.390 


-1-5 

-.2101 

-.1490 

7.9046 

-14.130 

1.02 

4 .8007 

-1.310 


-1.4 

-.2255 

-.1540 

7-7633 

-13.300 

1.04 

4 .7745 

-1.245 


-1.3 

- .2434 

-.1790 

7.6305 

-12.620 

1.06 

4.7496 

-1.200 


-1.2 

-.2622 

-.1880 

7-9041 

-11.9J0 

1.08 

4 .7296 

-1.140 


-1.1 

- .2834 

-.2120 

7 -3890 

-11.260 

1.10 

4 .7028 

-1.070 


-1 .0 

-.3064 

-.2300 

7.2724 

-10.790 

1.12 

4 .6814 

-I.060 


--95 

-.3188 

-.2480 

7-1649 

-10.210 

1.14 

4.6602 

-.990 


-.90 

-.3518 

-.2600 

7.0628 

-9-680 

1.16 

4.6404 

-•990 


-■85 

-.3453 

-.2700 

6.9660 

-9.220 

1.18 

4 .6214 

-.885 


-.80 

-■3596 

-.2860 

6.8730 

-8.810 

1 .20 

4.6037 

-.875 


-■75 

-.5742 

-.2920 

6 .7897 

-8.4io 

1 .22 

4 .9862 

-.845 


-.70 

-.3899 

- .3140 

6 .7016 

-8.030 

1.24 

4.9693 

-.800 


-.65 

- .4058 

-.3180 

6 .6213 

-7.610 

1.26 

4.5533 

--770 


-.60 

-.4226 

-.3560 

6.9492 

-7.290 

1.28 

4.5379 

-.740 


-■59 

-.4400 

-.3480 

6 .4723 

-7.000 

1.50 

4.5232 

-.710 


-.50 

-.4580 

-.3600 

6 .4023 

-6.76O 

1.32 

4 .5090 

-.670 


-.45 

-.4768 

-.3760 

6.3347 

-6.4oo 

1.34 

4 .4956 

-.660 


-.40 

-.4963 

-.3900 

6.2707 

-6.140 

1.36 

4 .4824 

-.615 


-■35 

- . 5161 

-.3960 

6.2093 

-9.910 

1.38 

4.4701 

-.605 


-■30 

- .5365 

-.4080 

6 . 1902 

-5.690 

1 .40 

4.4580 

-.600 


-.25 

-.5574 

-.4180 

6.0933 

-5.460 

1.42 

4.4460 

-.560 


- .20 

-.5787 

-.4260 

6.0307 

-5.330 

1.44 

4.4348 

-.540 


-■15 

■ - . 6004 

-.4540 

5-9864 

-5.070 

1.46 

4.4240 

-.520 


-.10 

- .6223 

-.4380 

9.9397 

-4.840 

1.48 

4.4136 

i -.490 


-.05 

-.6445 

-.4440 

9-8873 

-4.600 

1.50 

4.4038 

\ -.470 


0 

-.6667 

-.4440 

9.8409 

-4.500 

1.59 

4.3799 

! -.430 


.05 

- .6889 

-.4420 

9 .7999 

-4.350 

1.60 

4.3585 

j -.394 


.10 

-.7110 

-.4380 

9 -7920 

-4.190 

1.65 

4.3308 

-.370 


.15 

-.7329 

-.4340 

9 -7101 

-4.050 

1.70 

4.3203 

--334 


.20 

-.7946 

-.4260 

9 .6696 

-3-860 

1-79 

4 .3036 

-.316 


■25 

--7759 

- . 4180 

9.6310 

-3-810 

1.80 

4 .2078 

-.288 


■ 30 

-.7968 

-.4080 

9.9929 

-3.600 

1.89 

4.2734 

-.262 


■ 35 

-.8172 

- . 5960 

9.9969 

-3-530 

L.90 

4.2603 

- .254 


.40 

-.8370 

-.3900 

9.9216 

-3.420 

1.95 

4 .2476 

- .230 


.45 

-.8565 

-.3760 

9-4874 

-3.270 

2.0 

4 .2361 

-.212 


.50 

-.8753 

-.3600 

9.4947 

-3.170 

2.1 

4.2149 

-.177 


•55 

-.8933 

- .3480 

9 .4230 

-5.090 

2.2 

4.1972 

-.167 


.60 

-.9107 

-.3360 

9-3921 

-3.000 

2.3 

4 .1805 

-.136 


.65 

-.9275 

- .3180 

9.3621 

-2.880 

2.4 

4 .1669 

- .129 


■70 

-.9454 

- . 3140 

9-5533 

-2.810 

2-5 

4.1540 

-.114 


•75 

-.9591 

-.2920 

9.3092 

-2.690 

2.6 

4 .1426 

-.096 


.80 

-■9737 

-.2860 

9.2783 

-2.670 

2.7 

4 . 1330 

-■093 


.85 

- .9880 

-.2700 

9.2516 

-2.940 

2.8 

4.1237 

-.084 


-90 

-1.0015 

-.2600 

9.2262 

-2.480 

2.9 

4.1153 

-.071 


• 99 

-1.0145 

-.2480 

9-2014 

-2.390 

3.0 

4.1082 

-.065 


1.0 

-1.0269 

-.2300 

5.1179 

-2.350 

5.2 

4.0952 

-.054 


1.1 

-1.0499 

-.2120 

9.1940 

-2.260 

3-4 

4.0843 

-.044 


1.2 

-1 .0711 

-.1880 

9.1314 

-2.100 

3.6 

4.0755 

-.036 


1.3 

-1 .0899 

-.1790 

9.1096 

-2.130 

3.8 

4 .0684 

-.054 


1.4 

-I.IO78 

-.1940 



4.0 

4 .0616 



1.5 

1.6 

-I .1232 

-1.1377 

-.1490 

-.1310 











1-7 

-1.1908 

-.1190 






1.8 

-1.1627 

-.1110 






1-9 

-1.1738 

-.1010 






2.0 

-1.1839 

-.0890 






2.2 

-1.2017 

-.0745 






2.4 

-1.2166 

-.0625 






2.6 

-I.2291 

-.0945 






2.8 

-1.2400 

-.0470 






3-0 

-1.2494 

-.0400 






3.2 

-1.2574 

-.0345 






3-4 

-1.2643 

-.0305 






3-6 

-1.2704 

-.0290 






3-8 

-1.2762 

- .0230 






4.0 

-1.2808 





54 









55 












57 


TABLE VI.- WING-ROOT SHEAR, BENDING -MOMENT, 
AND PITCHING-MOMENT COMPONENTS 


q m p, 

lb 

lb/deg 

r > 

in-lb/deg 

M Y,a,r> 

in-lb/deg 

^n,r’ 

lfc/g 

in-lb/g 

M Y,n,r> 

in-lb/g 

Wing-root values derived from equations - 

f t^-deg 

(62) 

(65) 

(68) 

(63) 

(66) 

(69) 

0 

10 

20 

50 

40 

50 

60 

TO 

80 

90 

100 

0 

6,553 

12,052 

16,878 

21,244 

25,283 

29,076 

32,680 

36,134 

39,467 

42,699 

0 

1.893.500 
3,372,200 
4,588,400 

5.627.300 
6,540,600 
7,360,700 

8.110.500 

8.805.300 

9.456.300 
10,071,800 

0 

-1,307,900 

-2,334,000 

-3,183,600 

-3,915,700 

-4,565,900 

-5,157,200 

-5,704,800 

-6,219,600 

-6,709,200 

-7,179,400 

0 

1,226 

2,067 

2,664 

3,103 

3,428 

3,670 

3,850 

3,982 

4,076 

4,139 

0 

548,200 

942,300 

1.239.700 

1.472.200 

1.659.200 

1.812.700 
1,941,000 
2,049,850 
2,143,100 
2,223,600 

0 

-338,990 

-571,300 

-736,600 

-857,100 

-945,900 

-1,011,700 

-1,059,980 

-1,094,700 

-1,118,500 

-1,133,500 
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TABLE VII . - SAMPLE CALCULATIONS PER UNIT OF CHANGE IN LOAD FACTOR 
[w = 110,000 pounds; x cg = -211 Inches; = - 768.5 inchesj 


(a) Wing -root angle of attack, wing lcw*d, and tail load 


© 

© 

© 

© 

© 

© 

© 

© 


© 

qi V 

lb 

ft 2 -deg 

*tr,r> 

lb/deg 

%,a r. 
in -lb /deg 

©-? 

I'd , r > 
16/g 

^Y,n,r^ 

in-lb/g 

© + ^ 

/SPr © 

©© + © 
_ AL,. 

’ 

lb/g 

0 

1 *r[«f 

3 

0 

0 

0 

0 

0 

0 

0 

00 

(a) 

(a) 

10 

6,555 

-1,507,900 

4,851 

1,226 

-338,990 

59,112 

8.0627 

54,061 

1,878 

20 

12,052 

-2,534,000 

9,014 

2,067 

-571,300 

58,574 

^.2793 

53,642 

2,716 

30 

16,878 

-3,183,600 

12,754 

2,664 

-736,600 

58,192 

2.9992 

53,285 

5,450 

40 

21,21*1* 

-3,915,700 

16, 147 

5,105 

-857,100 

57,910 

2.3478 

52,980 

4,040 

50 

25,285 

-4,565,900 

19,540 

5,428 

-945,900 

57,700 

1.9493 

52,712 

4,576 

60 

29,076 

-5,157,200 

22,565 

5,670 

-1,011,700 

57,544 

I.6788 

52,484 

5,052 

70 

52,680 

-5,704,800 

25,255 

5,850 

-1,059,980 

57,427 

1.4820 

52,280 

5,440 

80 

56,151* 

-6,219,600 

28,059 

5,982 

-1,094,700 

57,540 

1.3317 

52,102 

5,796 

90 

59,467 

-6,709,200 

50,754 

4,076 

-1,118,500 

57,277 

1.2129 

51,945 

6,uo 

100 

42,699 

-7,179,400 

55,554 

1 

4,159 

i 

-1,133,500 

57,255 

I.H63 

51,804 

6,592 


(b) Pitching and bending mciSients 


© 

© 

© 

© 

© 

© 

© 

r' 

in -lb /deg 

M b,n,r' 

in-lb/g 

@© + © 
dn 

@® + © 
_ *v 

An 

* 

,p ® 
in. 

(D © . 

G 

II 

»? 

x 100 , 
c 

percent 

0 

0 

0 

0 

} - 206.6 

b 299.5 

0 

1 , 895,500 

548, 200 

-10,884,200 

15,814,900 

- 201,3 

292.5 

3.399 

3 , 372,200 

942,300 

-10,559,200 

15,373,000 

- 196.8 

286.6 

6.286 

4,588,400 

1,239,700 

-10,284,900 

15 , 001,200 

-193.0 

281.5 

8.723 

5,627,300 

1,472,200 

-10,050,400 

14,684,000 

-189.7 

277.2 

10.840 

6,540,600 

1 , 659,200 

-9,846,200 

14,408,800 

- 186.8 

273.3 

12.700 

7,360,700 

1 , 812,700 

- 9 , 669,600 

14,169,800 

-184.2 

270.0 

14.368 

8 , 110,500 

1,941,000 

-9,514,500 

13 , 960,800 

- 182.0 

267.0 

15.779 

8 , 805,300 

2,049,850 

-9,377,300 

13,775,900 

- 180.0 

264.4 

17.062 

9,456,300 

2,143,100 

- 9 , 256,100 

13 , 612,600 

- 178.2 

262.1 

18.217 

10 , 071,800 

2 , 223,600 

-9,147,900 

13,466,800 

- 176.6 

260.0 

19.243 


a Inde terminate value. 
^Rigid-wing value. 





Vortex station 
Downwash control po 



Figure 1.- Plan form of wing-nacelle combination used in sample problem. 









(b) Additional airload bending moment per unit of shear at center line 

Figure 3«- Spanwise variation of additional airload shear and bending 
moment for sample wing at two values of qm^. 




Dynamic pressure q, Ib/sq ft 

Figure 4.- Variation of qm^ with dynamic pressure for several Mach num 

bers for sample airplane. 
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